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We present a possible definition of a mobifity gap for a many-body quantum system, in analogy to 
definitions of dynamical localization for single particle systems. Using this definition, we construct 
"corrected" quasi-adiabatic continuation operators. We show that these operators have the same 
locality properties as the ordinary quasi-adiabatic continuation operators do in the case of a spectral 
gap, and that they approximate adiabatic evolution in the region with a mobility gap just as the 
ordinary operators do with a spectral gap. Further, under an appropriate definition of a unique 
ground state (equivalently, an absence of topological order as defined in the text), we show how to 
introduce virtual fluxes and prove bounds similar to those obtained on an energy for the effect of 
inserting 27r-flux. Armed with these results, we can directly carry over previous results proven in the 
case of a spectral gap. We present a statement of a higher-dimensional Lieb-Schultz-Mattis theorem 
for disordered systems (however, the lack of translational invariance presents us from proving the 
vanishing of the gap but rather only lets us prove a weaker statement that either the gap becomes 
superpolynomially small or the expectation value of the flux insertion operator varies in a particular 
way) ; we present a proof of decay of correlation functions; and we present a proof of Hall conductance 
quantization under very mild density-of-states assumptions defined later. We also generalize these 
definitions to the case of a "bulk mobility gap" , in the case of a system with boundaries, and present 
a proof of Hall conductance quantization on an annulus under appropriate assumptions. 

Further, we present a new "optimized" quasi-adiabatic continuation operator which simplifies 
previous estimates and tightens bounds in certain cases. This is presented in an appendix which can 
be read independently of the rest of the paper as it also improves estimates in the case of systems 
with a spectral gap. This filter function used decays in time at least as fast as C'(exp(— t")) for 
all a < 1, a class of decay that is called "subexponential" (a more precise and tighter description 
of what is possible is below). Using this function it is possible to tighten recent estimates of the 
Hall conductance quantization for gapped systems [8] to an error which also decays subexponentially 
(again, more precise descriptions are below), rather than just as an exponential of a power. 

INTRODUCTION 

Recently, the use of Lieb-Robinson[2Hl! bounds, combined with appropriately chosen filter functions has led to 
significant progress in proving results about quantum many-body systems. Examples include the higher-dimensional 
Lieb-Schultz-Mattis theoremjS], where this combination of techniques was introduced, decay of correlation functions 
in gapped systems [3HS1 IS], an area law for entanglement entropy for arbitrary one dimensional gapped systems |TD], a 
simpler proof of Goldstone's theorem for gapped Hamiltonians|13j. and, most recently, a proof of Hall conductance 
quantization for interacting electrons without averaging assumptions [8]. 

However, these results suffer from one major limitation: they require a spectral gap. However, in many cases 
we would prefer to require, instead, the weaker requirement of a mobility gap: a gap to propagating excitations. 
For non- interacting systems, the concept of localization has been around since Anderson's early workjy. Recently, 
interest has arisen in the possibility of localization in interesting systems. Interesting results include the possibility of 
a many-body localization transition |6], and theorems proving many-body localization for certain interacting systems, 
albeit in a special case that can be mapped to a non-interacting system[7,. In the single particle case, one can define 
localization in different ways. One way is in terms of the properties of the single particle eigenstates, while a different 
way is in terms of the dynamics [TT] . For a many-body system, the concept of single particle eigenstates no longer 
makes sense. Thus, we need to seek another definition. 

In this paper, we present a possible definition of the concept of a mobility gap. Further definitions are required 
to specify how the ground state should be distinguished from other states (definition |6] below) , and to generalize 
these concepts to open boundary conditions with gapless edge modes, as would appear in a Hall system. Using the 
appropriate definitions, we show how to generalize the concept of a quasi-adiabatic continuation operator [21 [H] to this 
kind of system, while preserving the needed locality properties. Given these results, we then are able to directly carry 
over many of the results previously shown using spectral gap. Under appropriate assumptions, we prove exponential 
decay of correlation functions, a version of a higher-dimensional Lieb-Schultz-Mattis theorem, and Hall conductance 
quantization. 

Our mobility gap definition describes an assumption that the propagation of low energy excitations is very slow. 
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This assumption is stronger than the usual Lieb-Robinson bound. Lieb- Robinson bounds hold for very general classes 
of Hamiltonians (essentially, any lattice Hamiltonian with short-range interactions and a bound on the interaction 
strength) . These bounds were introduced in [2] ■ In [9] , the idea of shifting certain terms in the equation of motion 
in a way that maintained the norm was introduced to show that these bounds hold in a way that is independent of 
the dimension of the Hilbert space on each site. In [3], a more general description on arbitrary lattices was given, 
albeit with dimension-dependent bounds, and in [3], the dimension-independent bounds were presented for arbitrary 
lattices and an extension to interactions decaying slower than exponential was given (this extension will be used in 
the appendix of this paper) . 

In this paper, we also present a more general definition of a quasi-adiabatic continuation operator, that contains the 
previous definitions. Further, we define an exact quasi-adiabatic continuation operator in the appendix with improved 
decay properties in time (exactly matching adiabatic evolution while also decaying in time as an exponential of a 
polynomial of the time). The use of this operator significantly simplifies the error estimates, in particular in cases 
where we need a Lieb-Robinson bound for quasi-adiabatic evolution. For example, it significantly tightens the error 
estimates in the recent proof of Hall conductance quantization. The appendix can be read separately. 

DEFINITIONS OF A MOBILITY GAP 

We consider lattice Hamiltonians of the following form: we assume that H is a. sum of terms 

Ho = Y.Hz, (1) 
z 

where each Hz is supported on set Z, and obeys the following. First, the diameter of every set Z is at most R. 
Second, 

sup,^||iJz|| < J, (2) 

Z3i 

where the supremum is over sites i. We let V denote the number of sites in the system. Thus, ||-ffo|| l£ JV- We refer 
to R as the "range" and J as the "interaction strength" . 

We use dist(-, •) for a metric on the lattice; we measure distances between pairs of sites, pairs of sets, or a site and 
a set using the same function. The distance between a pair of sets is defined to be the minimum over pairs of sites 
in the pair of sets of the distance between the sites, and similarly for the distance between a set and a site. We use 
diam(-) to indicate the diameter of a set. For any set A, we use A to denote the complement of A. We use L to 
denote various measures of the linear size of the system: for the case of a torus later, for example, we will consider 
an L-hy-L torus so that V = L^. 

We use to indicate the ground state of Hq, and similarly we use Pq = |^'o)(^ol to indicate the projector onto 
this ground state. We use Eq to denote the energy of state '^q. We use || • || to denote the operator norm and || -Hi 
to denote the trace norm and we use | ■ | to denote the P norm of a vector. 

We use C to refer to numeric constants of order unity. If we need multiple constants in the same expression, we 
use Ci, C2, .... We use poly(...) to refer to quantities bounded by a polynomial in their arguments. We use "computer 
science" big-0 notation: that is, indicating that a quantity is 0{x) indicates that it is bounded by a constant times x 
for sufficiently large x. We use exp(— poly(L)) to indicate that a quantity is C'(exp(— L"), for some a > 0. When we 
express bounds in term of the quantities L, Xmm, tmaxiT, this is always at fixed value of the quantities J/7, R, cioc, C 
(these quantities are defined below), and we use c to denote quantities which may depend on J/7, i?, c/oc, That is, if 
we state that a quantity is bounded by (J/Xmin) exp(— cL), we mean that the constant c is positive but may depend 
on J/^,R,cioct£,- When we state that a quantity is "superpolynomially small", we mean that it is superpolynomially 
small in L, for fixed J/7, R, cioc,^', in all such cases where we use the term "superpolynomially small" , we assume (and 
explicitly state) a polynomial dependence of quantities Xmin,T on L and a superpolynomially dependence oitmax on 
L. 

Before the definitions, some discussion is in order regarding "filter functions" . These functions play an essential 
role in the application of Lieb-Robinson bounds to many-body systems. The combination of these functions with 
Lieb-Robinson bounds was introduced in j9j. Broadly speaking, there are many places where, for a given operator O, 
we would like to construct a state /{Hq — Eq)0\^o), where /{Hq — Eq) is some function of the Hamiltonian Hq. That 
is, if Hq has eigenvectors with corresponding eigenvalues Ei, then /{Hq — Eq) has the same eigenvectors but has the 
eigenvalues /{Ei — Eo). In many such cases, the function /(w) that we would like to construct is not smooth near a; = 0. 
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The two functions that we would most hke to construct are the step function and the function l/w, which are used in 
proving correlation decay and in defining quasi-adiabatic continuation, respectively. If a system has an energy gap, then 
we can define a smooth filter function, /'(w), with the property that /'(w) is smooth and such that f{u)) — f'{(jj) is small 
for \lo\ larger than the energy gap. The smoothness property is used to show that the Fourier transform of /' is rapidly 
decaying in time, and hence to approximate /(i?o — £'o)0|^'o) = / dt/(t) exp(ii7oi)'^6xp(— «iJoi)|\E'o), by a local 
operator acting on VPo, using Lieb-Robinson bounds to show locality of exp{iHot)0 exp{—iHot) for fixed time. The 

smallness of /(w) — /'(w) for uj larger than the energy gap suffices to show that /{Hq — i?o)0|^o) ^ f{Ho — Eq)0\'^o) 
is small. There have been two main classes of filter functions considered. One can consider filter functions which 
decay exponentially in time at the cost of an exponentially small error in /(w) — /'(w). These functions, which we 
will call "Gaussian filters" (they are not equal to Gaussians, but have similar decay properties) often give the best 
bounds. The other class is filter functions was first considered by Osborne in as a modification of the Gaussian 
idea. We will call these functions "exact filters" . Exact filters have /(w) — /'(w) identically equal to zero for |a;| larger 
than the gap. These functions are easier to work with, but they often give bounds that decay only faster than any 
power (in the appendix, we present a construction of these filter functions that leads to the Fourier transform of / 
decaying as an exponential of a polynomial in time, for polynomial arbitrarily close to linear, but in the main text 
we content ourselves with superpolynomial decay); they also make it especially easy to prove Lieb-Robinson bounds 
for evolution under quasi-adiabatic continuation. In this paper, we will consider many of the definitions in generality, 
using abstract fiter functions. This will enable us to either find tighter bounds, or to simplify the proofs, depending 
on preference. 

Using such filter functions, one can define a quasi-adiabatic continuation operator to be an operator 

iV{H„dsH,) = J dtF{-ft)eM^Hst){dsH,)cxp{-iH,t), (3) 

where F{t) is some filter function such that its Fourier transform F{uj) approximates — l/w for |a;| > 1, F{t) decays 
rapidly in time, and F{t) is odd in time and F{0) = 0. Then, the Fourier transform of F['jt) approximates — l/oj for 
\uj\ > 7 and so we approximate adiabatic evolution given a spectral gap 7. Later, we modify this definition to account 
for a mobility gap. 

First, some definitions: 

Definition 1. For any set A, we define bi{A) to be the set of sites within distance I of set A. 

We use a Lieb-Robinson bound in the following form: 
Lemma 1. Given any operator O supported on a set A, for any I and any t with \t\ < 1/vlr, the operator 

0{t) = exp{iHQt)0 exp{~iHot) (4) 
can be approximated by an operator Oi (t) supported on bi (A) up to an error 

\\Oi{t)~Om<'^g{l)\A\\\0\\, (5) 
and also, for any operator U whose support does not intersect bi(A), we have 

||[O(i),t/]|l<^5(0HI|O||||C/|l, (6) 

where \A\ denotes the cardinality of the set A and g{l) decays faster than exponentially in l/R. The quantity v^ji 
depends on R, J, while the function g depends only on R. 

Proof. This is a minor variant of Lieb-Robinson bounds proven previously. See the appendix for an example of how 
such bounds are proven for a more general class of Hamiltonians. □ 

For short times, having the factor of \t\ in the above bound is useful, as it will help us deal with cases later that 
would otherwise lead to divergent integrals at short times. 

Definition 2. Given any operator O and function G, we define W^^g{0) to be the operator O filtered below 
energy 7 by 

VF^ ^(0) = 7 J dtGi-ft)exp{iHot)Oexp{-iHot), (7) 
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where the filter function G{t) and its corresponding Fourier transform G{uj) are chosen to have the properties that 
G{uj) is close to for > 1, G{uj) is close to 1 for \uj\ < 1/2, and G(t) is an even function of t and decays rapidly 
in t. 

Two specific examples of such, filter functions arc the following. First, we can define 

G{t) fq{t) ^ exp(-tV2g)[exp(i3t/4) - eyi^{-iZt/4)]/it. (8) 

For q — > oo, the Fourier transform of this is a filter onto frequencies between —3/4 and +3/4. That is, it is equal to 
unity for |a;| < 3/4 and equal to zero for |w| > 3/4. For finite g, fq{t) and its corresponding Fourier transform 
obey the following properties: 

• |/q(w)| < exp(— Cq) if \uj\ > 1, for some numeric constant C. 

• |/q(w) — 1| < cxp(— Cg) if < 1/2, for some numeric constant C. 

• ^ G exp{—t^ /2q), for some numeric constant G. 

This is an example of a Gaussian filter. Second, we can define an exact filer. An exact filter is a function G{uj) = 
Fiowit), where Fiow(t) is an even function, decaying faster than any power of t, with Fiowii^) = 1 for |a;| < 1/2 and 
Fiow{oj) = for |a;| = 1. 

If G{t) decays rapidly as a function of t, the operator W^^g(O) has the following localizability property which 
follows from a Lieb-Robinson bound: 

Lemma 2. If O is supported on set A, then for any I, the operator W^^g{0) can he approximated by an operator 
q{0) which is supported on bi{A) up to an error 

\\W,MO)-W\,am\<{ f \G{u)\du + g{l)\A\ fdu\G{u)\}\\0\\. (9) 



\y-\>h/vLR 



Proof. Let O' be defined by 



By a triangle inequality. 



dtG{-ft) exp{iHot)0 exp{-iHot). (10) 



10' -Oil < 7 / dt|G(7t)||10|' 

J\t\>l/VLR 



(11) 



Then, define Oi by 



By a triangle inequality and the Lieb-Robinson bound. 



\t\>l/vLR 

du|G(M)|||0|| 

\u\>lj/vLR 



Oi=l i ' dtG(^t)Oi{t). (12) 



|0,-0'|1 < 7 r'^'^ dt\G{^t)\g{l)\A\\\0\\ (13) 

-I/vlr 



< 7y dt|G(7t)|5(0l^ll|O|| 

< / du|G(u)|5(Z)|A|||0||. 



Eq. (|9| follows from Eqs. ( 11|13 | by a triangle inequality. □ 
As a corollary of the above result, we find that 
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Corollary 1. If O is supported on set A, then for any I, the operator W^j (O) can be approximated by an operator 
J (O) which is supported on bi{A) up to an error 

l|W^7./,(0)-</,(0)|| < {Cj^;^eM~ih/vLBf/2q]+C^gil)\A\)\\0\\. (14) 

and also 

Corollary 2. If O is supported on set A, then for any I, the operator Wj,Fia^{0) can be approximated by an operator 
W!^j^^^^(0) which is supported on bi{A) up to an error bounded by |yl||jO|| times a function decaying faster than any 
power of I. 

Definition 3. A Hamiltonian H is said to have a mobility gap 7 and localization length ^ and localization constant 
cioc up to time tmax if, for any operator O supported on set A and any filter function G, and any t with \t\ < tmax, 
there exists an operator W^°Q(0,t) with the following properties. First, for any I, VK^°Q(0,t) can be approximated by 
an operator supported on bi{A) up to an error in operator norm bounded by 

Qocexp(-//0||W^^,G(O)|| +niax|„l>^|GH|||0||. (15) 

Second, we require that the state produced by acting with W^^G(0)(t) = exp(iiJo0^7,G'(0) exp(— ii/o^) on the ground 
state is equal to the state produced by acting with the operator WlfQ{0,t) on the ground state, i.e., 

W-,.,{0)m^^W\%{0,t)^^. (16) 

Third, we have 

\\W\°;a{OM<\\W,,G{0)\\. (17) 

The above definition is our many-body version of tlie single particle definition of localization. It is an analogue of 



the definition of dynamical localization^^. The second term on the right-hand side of Eq. (15) reflects the "leakage" 
of states above the mobility gap due to the approximate nature of the filtering. 

It may happen that for a given system, there are several different choices of 7, ^, cjoc, tmax for which the system has 
a mobility gap. For example, in [7j, a very strong form of many-body localization was shown: for any operator O 
(even without filtering), the operator exp{illot)0 ejq){—iHot) could be approximated by an operator on a distance I 
with small error for I that was only logarithmically large in the time. Hence, by taking a given ^, for the system in 
[7], one can find a t^ax that is exponentially large in 

We assume that the ground state of Hq has energy Eq = 0. We define Xmin to be the second smallest eigenvalue of 
H. We will assume later only very modest requirements on Xmin- We will need tmaxXmin to be large, to control errors 
in quasi-adiabatic continuation. We will also need a unique bulk state as defined below: to have this, it suffices, but 
is not necessary to have Xmin > l/poly(i). This is a very weak requirement; for example, in a single particle system, 
the eigenvalue distribution is smooth, and so Xmin will typically be of order 1/V. 

We now define a corrected quasi-adiabatic continuation operator using the above definitions. 

Definition 4. Given a parameter- dependent Hamiltonian, Hq, an operator O, and functions F{t),G we define the 
corrected quasi-adiabatic continuation operator at mobility scale 7 and low energy cutoff to be the operator 
'D{Hq,0) defined by 

iV{Ha,0) ^ J F{2-ft)exp{iHot){0-W^,G{0))exp{^iHot)dt (18) 
+ [ F{X<t)Wi.%{0,t)dt, 



where the function F(t) has the property that its Fourier transform F(ijj) obeys 

F(0) = 0, (19) 



and where F is an odd function of time so that V is Hermitian. 
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Given a parameter dependent Hamiltonian Hs = Hz{s), we define 

V,=V{Hs,dsH,). (20) 

We also sometimes write T>f — 'D{Hs,dsHz{s)), so that 

2?,=^2?f. (21) 
z 

Definition 5. For each such quasi- adiabatic continuation, we define a function C{uj) by 

C{oj) ^ -Lf(u;/X<)G{uj) + i-F(^/27)(l - G(c.)). (22) 

The operator V wiU be used to approximate adiabatic evolution in a local way. To do this, we will require that 
C{io) be close to — l/cj for \uj\ > . Sec lemma ([6| and lemma ([t]) where we will show that 

9o|*o(s)) - iV^{Ho,dsHs)\'^o)\ < ||9.i/s||niax„|C(^) + 1/lo\. (23) 

The corrected quasi-adiabatic continuation operator here differs from previous ones, such as Eq. ([3|, by the addition 
of the terms involving W^°^(9s^s)- This will be used to account for low frequency components below the mobility 
gap. Our idea is as follows: for frequencies above the mobility gap, we use the large frequency to enable us to 
approximate the adiabatic evolution by an integral over a short-range of times, and hence with a local operator, while 
for frequencies below the mobility gap, we greatly increase the time scale we use to approximate adiabatic evolution, 
but we use the assumption of localization below the mobility gap to keep the operators local. 

One example of such a function F gives a Gaussian corrected quasi-adiabatic continuation operator. In [5], the 
quasi-adiabatic evolution operator was defined by 



— du dtexp{iHot)Oexp{-iHot)cxp{-u'^/2a'^), 
av27r J Jo 



for some parameter a, while in [9j I12j a more complicated integral was used. Eq. (24) can be re- written as 



(24) 



J dtF{t) exp{iHot)0 exp{~iHot), (25) 

where F{t) is defined to be 

i r°° 

duexp(-uV2a^)> (26) 



av 27r 

for t > and F(t) is an odd function. We can use this function F(t) in our corrected quasi-adiabatic continuation 
definition, using G = fq for the filter function, getting 

|C(w) + < C(l/A™„) exp[-C(A</A™„)2a2] + C(l/A„„„) exp[-C<z] + ^(1/7) exp[-CaV2], (27) 



for |w| > Xjuim a-nd we have 



\F{t)\<Cexp[-t'^/2a^]. (28) 

Alternatively, we can define an exact corrected quasi-adiabatic continuation operator at mobility scale 7, and low- 
energy cutoff A< by a function F{t) where F{t) is some function which has the property that its Fourier transform, 
F{ijj), is odd and infinitely differentiable and equals —l/co for \uj\ > 1. Then, we have 

C{u)^-l/uj (29) 

for |a;| > A^. Since F{uj) is infinitely differentiable, F{t) decays faster than any power of t. 
Finally, we need one more definition ^23]: 
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Definition 6. We say that a Hamiltonian H has an (l, t) unique ground state if the following holds for all e > 0. 
Given any density matrix p such that, for all sets A with diam(yl) < I the inequality 

|iTr^(p-Po)||i <e (30) 

holds, then 

Hp - Polli <TVi. (31) 

We now show that give a bound on the smallest eigenvalue, \min, then an (/,t) unique ground state follows for a 
T that depends on Xmin- 

Lemma 3. // the second smallest eigenvalue is at least Xmin for a Hamiltonian H with range R and interaction 
strength J, then it has an (R, 1\J JV/Xmin) unique ground state according to the above definition. 



Proof. Suppose Eq. (30) holds. Note that Tr(pi7) = J^z^^'^iP^z)- Then, since each Hz is supported on a set of 
diameter R, we have that Tr((p - Po)Hz) < e\\Hz\\. Summing over Z, Ti{pH) - Eq) < \\H\\ < eJV. We now 
maximize the 1-norm difference between p and Pq subject to this constraint on the energy difference. The maximum 
is obtained when p = \tp){'ip\, for V' — cos((?)5'o + sin(6')^'i, with "ifi being an eigenstate of H with energy Xmin- Using 
the estimate of tr{pH), sin(0)^ < eJV/ Xmin- In this two dimensional subspace, p — Pq equals 



-sin(6i)2 cos(6') sin(e') 
^cos(6') sin(6l) sin(6l)2 

and Hp- Poll 1 = 2^sin(6i)4 + cos(6')2 sin(6l)2 = 2 sin(6') < 2^ejy/A„„„ 



(32) 



Thus, Eq. (31 1 follows with t = 2^/jV/X~- □ 



One fundamental idea in and [8], was to show that quasi-adiabatic evolution around certain closed paths in 
parameter space left the energy almost unchanged at the end of the path; then, using the existence of a spectral gap, 
the fact that the energy was almost unchanged was used to show that we had returned to a state close to the ground 
state. Here, we will use this {1,t) unique bulk state assumption instead of a spectral gap since we may not have a 
spectral gap. In fact, however, since we will only need the {1,t) assumption for r > l/poly(L), from lemma ([3| it 
suffices to have a minimum eigenvalue Xmin which is greater than or equal to l/poly(i). So, the (/,t) unique ground 
state assumption follows from a very weak assumption on Xmin, as claimed above. 

This unique ground state assumption is physically necessary when we prove Hall conductance quantization later 
under a weaker assumption of a mobility gap, compared to the spectral gap assumption of [S]. On physical grounds, 
we need to have the {1,t) unique ground state for the following reason: consider a fractional Hall system on a torus 
with multiply degenerate ground state and then a spectral gap (not just a mobility gap) to the rest of the spectrum. 
This system will not display integer Hall conductance quantization. However, it will actually show a mobility gap, 
up to exponentially large times t, since after filtering any operator acting on one ground state can only have matrix 
elements to one of the other ground states, and since the splitting between the ground states is exponentially small 
in system size, the filtered operator will be almost unchanging in time. Thus, we do not expect that there is any Hall 
conductance quantization theorem in the absence of some condition like the (/,t) unique ground state condition. 

As further justification for our [I, r) unique ground state definition, we note that if this definition does not hold for 
I — L/2 — 1 and r << 1, then there exists another state ip orthogonal to with the property that, given any local 
operator O with support on a set of diameter less than half the system size, if the operator O is projected into the 
two-dimensional spanned by \E'o,V' it is close to the identity operator. This is a definition of topological order (see, 
for example the definition of (Z,e) topological order in [H]). Thus, since we will only use this unique ground state 
definition for I ~ const, x L later, we are in fact only requiring the absence of topological order. 

A further reason for introducing the unique ground state assumption is that later in the context of the Hall effect 



with boundaries we will need a different unique bulk state assumption, definition (10), which generalizes this unique 
ground state assumption. 



CORRELATION DECAY 



The most basic result to show using these definitions is the exponential decay of correlations in a system with an 
unique ground state and a mobility gap. We do all these calculations with Gaussian filter functions. This section can 
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be read separately from the later sections of the text, because it only relies on the assumption of a mobility gap and 
does not use the definitions of quasi-adiabatic continuation operators. 
We begin with a lemma: 

Lemma 4. Let -P-y/2 denote the projector onto eigenstates with energy greater than or equal to Eq +7/2. Let Oa, Ob 
be operators supported on sets A, B with dist(A, i?) = I. Suppose that Ob has that property that 



(33) 



(1 - P^/2)Ob*0 
for some S. Then, 

\{^o,OaOb^o) - {^o,Oa^o){-^o,Ob'^o)\ < c[exp{^Ch/2vLR) + mini\Al\B\)g{l) + s}\\Oa\\\\Ob\\. (34) 



Proof. The proof basically follows previously proven correlation bounds [3] . Assume without loss of generality that 
(5'o, Oyi^o) — {'^0: Ob'^o) — 0- For any operator X, following 5,, we define X+ by 



hm — / diexp(ii7ot)Xexp(-zi/ot)^^^^^L(^^^^/^l!M 
e-i.o+ 27r / it + e 



(35) 



for some q which will be chosen equal to lj/2vi^j^ below. 

The operator is equal to X{t) — exp{iHot) convolved against the function (l/27r)(exp[— (t7/2)^/2g])/(it + e). 
For q 00, the Fourier function converges to a step function, vanishing for negative oj, and unity for positive w. For 
finite q, one may show that for any such operator X, we have 



as shown in [5], and also that 



X-X+]R 



7/2 



/27rq 



^^ cxp(-g/2) 



/27rq 



Further, one may show from Eq. ( 33 I that 

(1-P^/2)(Ob|*o)-0+|vI/o) 

Thus, 



< 6. 



(6+ -Ob) 1*0) 



(36) 



(37) 



(38) 



(39) 



We now estimate the commutator [Og , Oa] by the Lieb- Robinson bound and the usual trick of splitting the time 
integral into integrals over early times (|t| < 1/vlr) and the integral over late times {\t\ > 1/vlr)- By a triangle 
inequality 



m,OA]\\ < ^/d^ ^-p'-y/^^% o.(^),o.]ii 



(40) 



< 



< 



< 



1 

2n 



2tt 



\t\>l/vLR 



\t\>l/vLR 



d. °''-y''^'' ii[o.(.),o.iiKi 



\t\<l/vLR 



\t\<l/vLR 

c{iqvlJ{f^'j)eM~il'lV^<H)/H + \B\g{l)}\\OA\\\\OB\\. 



dt'^\B\gm[OBit),OA]\\ 
dt'^\B\gil)\\[OBit),OA]\\ 



Using the triangle inequality | (*o, OaOb^-q) | < K^-o, Oa{Ob - 0+)*o) | + I (*o, 0+Oa*o) I + II [Oa, 0+] \\ , choosing 
q ^ l"f /2vlr, '^c find that 



\{'fo,OAOB'fo)\<c[^exp{-Ch/2vLR)+gil)\B\+6}\\OA\\\\Oi 



(41) 
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Since (vlr/Ij) exp{-Clj/2vLR) < C exp{-lj/2vLR) for I > vlr/j, we have 

\{^o,OaOb'^o)\ < c[eM'Ch/2vLR) + g{l)\B\ + 6}\\0a\\\\0b\\. (42) 
One can replace the \B\ with |^| by applying the Lieb- Robinson bound to the evolution of Oa rather than Ob- 



□ 



max • 
min • 



Theorem 1. Assume that H has a mobility gap 7, localization length ^ and localization constant cioc up to time t 
Let Oa, Ob be operators supported on sets A, B with dist(^, B) — I. Assume that the lowest eigenvalue of H is A 
Then, 

|(*0OaOb*0)| < [Cg{l/2)\A\+eM-tmaxXrn^n)+l0gi2vLRtrnax/l)[cloceM-l/0+C<iM-Ch/^VLR))^ 

(43) 

Proof. Assume, without loss of generality, that (Oa) = {Ob) = 0. Define Wjj^{Ob) to be the operator Ob filtered 
below energy 7 as before. Define Z = W^jJO), with q = Ij/2vlr. Then, 

(*o, OaOb^o) = (*o, Oa (Ob - ^) *o) + ('J'o, Oa (z - W^j^ (O)) ^o) + (*o, OaW^^,/, (0)*o) • (44) 
Let X = 0b- Z. We have 

(l-P^/2)|XvI/o)| < \\wlll{OB)-W^.j^{OB)\\+eM-Cq)\\OB\\ (45) 
- i^ h/2vLR ''M-illl'^VLRf/H + Cg{l/2)\A\ + exp(-C<7)) ||0b||. 
Further, X is supported on a set which is at least distance 1/2 from A. So, by the previous lemma, 

\{-^o,OaX^^)\ < c[exp{~Ch/2vLR) + \A\g{l) + Cg{l/2)\A\+exp{-Cq))}\\OA\\\\OB\\ (46) 
< C[exp{-Cq) + \A\g{l/2))}\\OA\\\\OB\\. 

Note that if desired, the term |A| in the above expression can be replaced by the cardinality of the set of sites within 
distance 1/2 oi B and the bound still holds. 

We now consider the term (^o, Oa(^Z - W-yj^(0)^*o) in Eq. (44 1. This is bounded by 

\\w'^{1{Ob)-W^j^{Ob)\\ (47) 



< exp[-{l^/2vLR)y2q] + Cgil/2)\A\) ||0, 



l-i/2vLR 

We finally consider the term (^'o, OaW^7./5 (0)5*0). Define Y by 



y = / <^t\t\<t^.w\:uo,t)d{t) (48) 

' j dt|t|>t„„, exp(ii?oOW^7 J, (O) exp{-iHot)d{t), 



where the function d{t) is defined by 



d{t) = exp[-(tA^,„)V2gi] '''P^'^.T^^ \ (49) 

it 



for some qi. This is the similar to (35), but we use the filter d(t). For (ji — >■ oo, the filter d{t) has a Fourier transform 
equal to unity for < w < 2jt and vanishing for w < or a; > 2'yt and so at infinite qi, \Wjj^{0)'^() — F^o| < 
|P27W^7j^(0)^o| ^ exp(— Og)||0||. For finite qi, the filter d{t) approximates this filter. The filter d{t) is chosen not 
to have a singularity at t = 0. So, using the assumption on Xmin, we find that 

W^7./,(0)|*o) - i^|*o)| < (exp(-Og) + "'^'^g/^^ )||0b||. (50) 
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Similarly, 



[^o\Y < (exp(-C,) + '"pJ^/^V bII. (51) 



We now estimate the commutator ||[F, 0^]||. We do this using a triangle inequality. The integral over \t\ > t^ax in 
(|48| is bounded by C{qi/tmax>^min) exp[-(t™aa;Amm)^/2qi]. We break the integral up of |i| < tmax into two different 
parts. First, a part with \t\ < 1/2vlr- Second, a part with l/2vLfi < \t\ < tmax- 

We now bound the integral over \t\ < 1/2vlji. By the localization assumption, 

\\[W,j^{OB,t),OA]\\ < [ciocexp{-^l/0 + Cexp{-Cq))\\OA\\\\OB\\. (52) 
For |t| < l/2vLjf , d{t) is bounded by 27. Thus, the integral over \t\ < 1/2vlji is bounded by 

hl/vLR) (qoc exp(-VO + Cexp(-cg)) ||0a|| ||Ob|| . (53) 
We next bound the integral over 1/2vlj^ < \t\ < tmax- By the localization assumption, 

\\[W,j^{OB,t),OA]\\ < (qo, exp(- + C cxp(-Cg)) II II II Oi3 II. (54) 
Thus, the integral over 1/2vlj^ < \t\ < tmax is bounded by 

log{2vLRtmax/l)(^ciocexp{~l/^) + C exp{-cq)'^\\OA\\\\OB\\. (55) 
We pick qi — tmaxKnin- Thus by a triangle inequality, 
\{^o,OaZ^o)\ < {Cgil/2)\A\+exp{ )+\og{2vLRtmaxll)(^ciocex-p{-l/C)+C ex^i-Cl-f I2vlr)^^\\0a\ 

Thus, 

|(*0OaOb1'0)| < \Cg{l/2)\A\+exp{-tmax\nzn)+^Og{2vLRtmaxll)(^CloceyiV{-l/0+C 



ObI 
(56) 



Os||. 
(57) 

□ 

PROPERTIES OF CORRECTED QUASI-ADIABATIC CONTINUATION OPERATOR 

We now consider the properties of the corrected quasi-adiabatic continuation operator. There are three basic 
properties used previously in studying these systems. First, the quasi-adiabatic continuation operator should be local, 
in that it should be a sum of operators, each of which is exponentially decaying in space, in the sense that each 
such operator can be approximated to exponentially good accuracy by an operator with finite range. Second, the 
quasi-adiabatic continuation operator should approximate the exact adiabatic evolution of a state in a region in which 
the system has a mobility gap and a sufficiently large Xmin- Third, the quasi-adiabatic continuation operator should 
produce the correct Berry phase: that is, we should have that 

(*o(s),2?(i/„5,i7,)*o(s)) =0, (58) 



where 5*0(5) is the ground state of Hg. This last property follows immediately from Eq. (191 in the definition of the 
corrected quasi-adiabatic continuation operator and was emphasized in 13] . We now show the other two properties 
for Gaussian and exact filter functions. 

Lemma 5. Consider a corrected quasi-adiabatic continuation operator and a Hamiltonian with a mobility gap. Define 
a function B(t) to be the convolution of F{jt) with 6(t) — ^G{'yt) (here, S{t) denotes the Dirac delta- function) . Then, 
for any operator O with support on set A and any U with support on set B with dist(yl, B) > I, then 

||[P(i/o,0),C/]||<i?(0||0||||C/||, (59) 
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where E{1) is defined by 

E{1) = {maxt{\Bm{l/vLR)\0\g{l)+ [ dt\B{t)\} (60) 

Qoc exp(— ;/^) +max|i^|>^ |GH|)+ / dt\F{X<t)\}. 



\t\>l/vLB. 

+ 

l\t\>t 

and further, for any I, there exists an operator O' with support on bi{A) such that 



\\V{H^,0)-0'\\ (61) 



Proof. We consider both terms in the definition of V, Eq. (18 1, separately. For the first term, we wish to bound 



We note that the operator 



is equal to 



\ j dtF{t)eM^Hot){0-W^,G{0))eM-iHot)M\V (62) 
j dtF{t) exp{iHot){0 - W^,g{0)) exp{-iHot) (63) 
J dtB{t) exp{iHot)0 cxp{-iHot), (64) 



for a function B{t) equal to the convolution of F{t) with S{t) — 76(7^). 
We use a triangle inequality: 



\\[J dtB{t)exp{iHat)Oexp{-iHot),U]\\ (65) 
< / dtB{t)\\[exp{iHot)Oexp{-iHot),U]\\ 
+ j dtB{t)\\[exp{iHot)Oexp{-iHot),U]\\. 

J\t\>l/VLR 

The first term is bounded using the Lieb-Robinson bound by 

mi,xti\Bitmi/vLR)\0\gil)\\0\\\\U\\. (66) 



The second term is bounded by 



We now consider the term 



/ dt\BmO\\\\U\\. (67) 

J\t\>l/VLB. 



\\[J dtF{X<t)w!^%{0,t),U]\\ (68) 

< / dtF{X<t)\\[Wl;,^aiO,t),U]\\ 
J\t\<t„,„„. 



i\t\<t„ 

+ f dtF{X<t)\\[Wi^%{0,t),U]\\. 
To bound the integral over \t\ < tmax, we use the localization assumption, to bound this by 



tniax{ciocexp{-l I S,) + max|„|>.^ |G(w) |^ ||0||||C/||. (69) 
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The integral over \t\ > tmax is bounded by 

dt|^^(A<<)|||0||||[/||. (70) 

|t|>t„ax 

Putting these results ( 66|67|69|70 ) together, Eq. (59) follows. Since Eq. (59) holds for all operators U, Eq. (61) 
follows: to see this, define O' to be 

O' = j AU UV{Hq,0)U\ (71) 
where the integral is over all unitary rotations over sites not in hi{A) with the Haar measure. This trick was introduced 

in m]. □ 

Applying these results to the Gaussian corrected quasi-adiabatic continuation operator and Gaussian filter function, 
we have the corollary: 

Corollary 3. For a Gaussian corrected quasi-adiabatic continuation operator and a Hamiltonian with a mobility gap, 
the error term E{1) in is bounded by: 



\mH^,0),U]\\ (72) 
< {c{l/vLB)\0\g{l) + C{h/vLB){l/vLR){o?eM~C{h/vLRf/a^]+qeM~C{l^/vLR)^/q]^^ 

max /«)V2]} 



x||0||||f/|| 



Proof. Note that F{t) decays as exp[— t/a^], while G{t) decays as exp[— Ct^/q], so B{t) is bounded by 
(74) ^exp[— C(7i)^/Q;^] + exp[—C (7^)^/9]^ Note that the function B{t) decays as exp[— C(i/a)^] for large t for some 
C. The rest follows immediately from the definitions. □ 

Similarly, for exact corrected quasi-adiabatic continuation operators we have the corollary: 

Corollary 4. For an exact corrected quasi-adiabatic continuation operator and a Hamiltonian with a mobility gap, 
the error term E[l) in l6(f\j is bounded by: 



E{1) < C{l/vLB.)\0\g{l) + -Qi{Ij/vlb) + Cloct^ax exp(-Z/C) + ^Q2[\'^t„,ax). (73) 

7 

where the Junctions Qi, Q2 decays faster than any power of their arguments. 

This implies the following "superpolynomial localizability" property: 

Corollary 5. The exact corrected quasi-adiabatic continuation operator 'D{H a, dsHz{s)) can be approximated by an 
operator supported on bi{Z) up to an error bounded by \Z\\\dsHz{s)\\ times I/7 times a function decaying superpoly- 
nomially in I plus 1/A^ times a function decaying superpolynomially in X^tmax- 

We now show that this corrected quasi-adiabatic continuation operator approximates the adiabatic evolution of 
states for Hamiltonians with a mobility gap. 

Lemma 6. For a Gaussian corrected quasi-adiabatic continuation operator T), and a Hamiltonian Hq with a mobility 
gap, and a given Xmin > 0, we have 

ds\^o{s))-tV,{Ho,dsHs)\^o)\ (74) 

< \F{u:/X<)fiu;,q) + F{ujh){l - f{io,q)) + lM\\dM 

< C((l/A™„)exp[-C(A</A™„)2a2] + (l/A„„) exph^g] + (I/7) expi-CaVS], ) P.i/,!]. 
where ^'o(s) is the ground state eigenvector of Hs and the partial derivatives are taken at s — 0. 
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Proof. Let ^'i(s) denote a complete basis of eigenstates of Hg, with corresponding eigenvalues Ei{s). We have 

dsMs) = E go(0)-g,(0) ^'^'^)<^-W, (a.i/.)vl'o(0)), (75) 

by linear perturbation theory. 

Also, by the definition of the corrected quasi-adiabatic continuation operator, 

iV^{Ho,d,Hs)'^o - ;EC(£;, -£;o)^'o(0)(^'^(0), (a,i/,)*o(0)), (76) 



ao|*o(s))-iI?a(Ho,5.i?,)|*o) < |19.H,|lmax„|C(c^) + l/t^|. (77) 



so by Eq. (27), Eq. (74) follows. □ 



The exact version of the above lemma is much simpler: 

Lemma 7. For an exact corrected quasi-adiabatic continuation operator V and exact filter function Fiow, and a 
Hamiltonian Hq with a mobility gap, and a given Xmin > we have 

d,\^ois)) ^iVs{Ho,d,H,)\^o)- (78) 



Proof. This proof is immediate from Eq. ( 29 ) . □ 



In some cases, the Gaussian operators above allow tighter estimates. However, from now on, for simplicity of 
estimates, and for the simplicity of expressing the results, we will use the exact quasi-adiabatic operators (further, 
using the construction in the appendix, one can see that in fact our bounds here, which will be expressed only as 
superpolynomial decay in L, in fact become "subexponential" decay in L, as defined in the appendix). In particular, 
the Gaussian operators give slightly better bounds in the Lieb-Schultz-Mattis case of the next section (we omit the 
results for simplicity), while the exact operators will actually lead to better bounds in the Hall conductance section. 

LIEB-SCHULTZ-MATTIS-TYPE THEOREMS 

We now consider applying these results to prove Lieb-Schultz-Mattis-type theorems. The ideas here will be needed 
for the Hall discussion later. We consider a system of linear size L. For definiteness, we consider hypercubic geometry; 
that is, we consider square geometry in two dimensions {L by L), cubic geometry in three dimensions, and so on. 

In this section, we will assume periodic boundary conditions in one direction of the hypercube, which we call the 
X direction. This does not mean that we assume translation invariance in that direction as done in [9l I14j . Rather it 
means that our metric dist(-, •) only measures distances between points mod(i). 

The one dimensional Lieb-Schultz-Mattis theorem, as later generalized in ^l8^, was a statement about translationally 
invariant one-dimensional quantum systems, with finite range and finite strength interaction, and with a conserved 
local charge. Having a conserved local charge means that there is some operator qi, defined on each site i, such that 
Q = qi commutes with the Hamiltonian, and such that qi has integer eigenvalues with \\qi\\ < qmax for some given 
Qmax- Then, assuming that Q is not an integer multiple of L, it was proven that the gap from the ground state to the 
first excited state decays as 1/L. 

This was generalized to higher dimensions in [S]. The most general statement, in |15) . is that if Q is not an integer 
multiple of L, then the gap from the ground state to the first excited state decays is bounded by 0{\og{L) / L). 

It is important to understand what being an integer multiple of L means. The work [9l [15^, only required trans- 
lational invariance in one direction, the x direction. However, if the system is an L-by-L' torus, and has a filling 
fraction Q/V = 1/2, then Q/L is non-integer if L' is odd. The restriction to odd width arises because we use ideas of 
flux insertion to construct a state which has low energy and which has a different momentum compared to the ground 
state, thus proving bounds on the energy gap variationally. The major improvement compared to the one-dimensional 
result was the ability to handle systems whose aspect ratio was of order unity. 

In this section we consider disordered systems, without translation invariance. Thus, we will certainly not be able 
to prove the existence of low energy excitations in this section, because in the absence of translation invariance there 
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exist Hamiltonians in which Q is not an integer multiple of L but with a unique ground state and a spectral gap. 
Instead, what we will show is the following. We will construct a flux insertion operator which inserts 27r flux in a 
vertical line, and apply it to the ground state. We will show that the expectation value of the energy of the resulting 
state is exponentially small in L. Thus, cither A,„i„ is exponentially small in L, or the flux insertion operator acting 
on the ground state produces a state which is superpolynomially close to the ground state multiplied by a phase. We 
then show in the latter case, where the flux insertion operator is superpolynomially close to acting on the ground 
state by a phase, that if Q/L is non-integer, then this phase depends in a particular way described below on which 
line is chosen for the flux insertion. 

To define the flux insertion operator, we need to define the Hamiltonian with twisted boundary conditions. Let Qx 
be defined by 



l<x(i)<L/2 

E 



Qx= V ft, (79) 



where x{i) is the x-coordinate of site i. That is, Qx is the total charge in the half of the system to the left of the 
vertical line with x = I//2 + 1 and to the right of a; = 0. Let 

H{ei,92) = Y,Hz{9ue2), (80) 
z 

where Hz{0i,92) is defined as follows. If the set Z is within distance R of the vertical line a; = 0, then Hz{0i,O2) = 
ex]){i9iQx)Hz exp{—i9iQx); if the set Z is within distance R of the vertical line x = L/2, then Hz{9i,62) = 
exp{—i02Qx)Hzexp{i92Qx)', otherwise, Hz{9i,92) = Hz- Note that, 

H{9,-9) = exi,{i9Qx)Hcxp{~i9Qx)- (81) 

We define an operator 2?,, to be an exact corrected quasi-adiabatic continuation operator with = Xmin describing 
quasi-adiabatic continuation along the path from 9i = 3,02 = —s. Wc have that dsH{s) = dsHz{s), and dsHz{s) 
is nonvanishing if Z is within distance R of the line at a; = or if Z is within distance R of the line at a; = L/2. 
Let ©(^'(s) denote the sum of terms in dsH{s) near the line at .t = and let O'^^'(s) denote the sum of terms in 
dsH{s) near the line at a; = L/2, so that = V{Hs,0'^^\s)) + V(Hs,0^'^\s)). Note that the Hamiltonians 
are all unitarily equivalent. So, if there is a mobility gap at s = 0, then there is a mobility gap for all s. Under 
the assumption of a mobility gap, because of the superpolynomial localizability property, we can approximate the 
operators 'D{Hg,0''^^) and 'D{Hg,0^'^^) by operators vi^^ and vi'^^ supported within distance less than L/8 of the 
respective lines a; = and x = L/2 up to superpolynomially small error. Note that the supports of D^^'^^ do not 
overlap. 

We define the fiux insertion operator, Wi, as follows. We define a unitary ui^^ by 

= S' exp{i /' ds'©^'^}, (82) 
Jo 

where S' denotes that the integral is s'-ordered. Then, we set Wi = U2^K We define W2 similarly: we define a unitary 
by 

[7,(2) = S' exp{i r ds'P^'^}, (83) 
^0 

(2) 

and we set W2 = ■ 
We define 

W = WxW2. (84) 

Now, we claim that: 

Lemma 8. Assume that the system has an {L/8,t) unique bulk state, with t greater than or equal to l/poly(i). 
Assume that the system has a mobility gap, with tmax superpolynomially large in L. Then, 

min^,,l^,|=i(|M/i|*o) - 2i|*o)|) (85) 

is bounded by J/7 times a function decaying superpolynomially in L plus poly(Z/)J/A^ times a function decaying 
superpolynomially in X^tmax- 
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Proof. Consider any operator O supported on a set of size A at most L/8. If A is not within distance L/8 of the hne 
X = 0, then 

{^q\wIoWi\^o) = (*o|0|*o), (86) 

since Wi commutes with O and Wi is unitary. 

On the other hand, if O is within distance L/8 of the hne a; = 0, then 

(*o|VF/01^i|*o) = {'i'o\W^WlOWiW2\'^o) (87) 
= {^o\W^OW\-^o), 



since W2 commutes with Wi and O. 
Define 



U = S'exp{i / ds'V,,}, 
Jo 



along the path s — 9i — —62- Then, the operator norm difference, \\U — W\\ is bounded by J/7 times a function 
decaying superpolynomially in L plus poly(L) J/A^ times a function decaying superpolynomially in X^tmax- 

However, U is an exact quasi-adiabatic evolution, and the Hamiltonians Hs are unitarily equivalent by Eq. (81 1. 
Thus, 

C/|*o) =^|*o), (89) 

for some z with \z\ = 1. Thus, | (5'o|W^^OVr|4'o) ~ (^'o|0|^o)| is bounded by is bounded by J/7 times a function 
decaying superpolynomially in L plus poly(L) J/A^ times a function decaying superpolynomially in X^tmax- 
Thus, from Eqs. ( 86|87 ), for any operator O, 

|(*o|l^/OM^i|^'o)-(*o|0|*o)| (90) 

is bounded by J/7 times a function decaying superpolynomially in L plus poly(X)J/A^ times a function decaying 
superpolynomially in X^tmax- 

However, the lemma then follows by the assumption of {L/8,t) unique bulk state. □ 

We need an estimate of the Berry phase. This is similar to the ideas in[T2]. Above, we claimed that 

C/|*o) =^|*o), (91) 
for some z with \z\ = 1. We now determine the value of z. We claim that 
Lemma 9. 

z ^ exp{-i2TrQx), (92) 

where 

Qx = (*o|Qx|*o). (93) 



Proof. By Eq. 

i^slVsl'^s) =0- (94) 
However, since Vs an exact quasi-adiabatic evolution operator 

= i(Ox-c)*„ (95) 



for some constant c. Eq. (94| lets us determine c so that Wsl'^s) = (Qx — Qx)^s- Thus, U\'^o) = 
exp{i2TrQx) exp{—i2TTQx)^o = exp{—i2TTQx)'^o- n 
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We now consider the phase zi in lemma Consider a flux insertion operator Wi{xo) defined precisely as the 
above Wi was defined above, except for inserting the flux along the line x = xq, rather than along the line a: = 0. 
That is, we define 



l+xo<x{i)<L/2 

E 



Qx{xo)= V q^, (96) 



and 

Hiei,e2,xo) = Y.Hzi9i,e2,xo), (97) 

z 

where Hz{di, 02, xq) is defined as follows. If the set Z is within distance R of the vertical line x = xq, then Hz{di, O2) = 
exp(i6'i(3x)-ffz exp(— i6'2Qjs:); if the set Z is within distance R of the vertical line x — L/2, then Hz{0i,02) = 
exp{i92Qx)Hz exp{—i62Qx), otherwise, Hz{Oi,02) — Hz- Then, we can define v''s'\xq) and Wi{xq) simflarly to 
before. Note that M^i(O) = Wi. 



For any define 



Q^= (*o|g.|^o). (98) 

Then, we claim that 
Lemma 10. The quantity 

|(*o|VFi(0)^VFi(l)|*o) - cxp(-i2^Qi)l (99) 

is hounded by J/7 times a function decaying superpolynomially in L plus poly(L)J/A^ times a function decaying 
superpolynomially in X^t^ax- 

Proof. We have 

(*o|W^i(0)^W^i(l)|*o) = (*o|W^i(0)V2V2M^i(l)|*o). (100) 

By the above lemma, W2W^i(0)|\I'o) is superpolynomially close to '^0 times some phase ^(0). Similarly, one can 
prove M^2W^i(l)|^o) is superpolynomially close to ^Pq times some phase ^(1); note that the support of the operators 
Wi(l) and W2 still do not overlap. Thus, 

(*o|Tyi(0)^Tyi(l)|*o) (101) 

is superpolynomially close to 2;(0)^z(l). 

However, by the above lemma, z(0) is superpolynomially close to exp(— z27r ^^^-^ Q^) and z(l) is superpolynomially 
close to exp(— i27r^^£2 Qx)- So, zi(0)'''2;i(l) is superpolynomially close to exp{—i2T:Qi). □ 

Note that the choice of the line a; = was arbitrary. We can pick any line x = xq to insert fiux into, and in that 
way define an operator Wi{xo), and then we can choose another line x = xq + L/2 to define W2{xq) and use those 
two lines and repeat the above proofs, showing that the quantity 

\{^o\Wi{xo)^W,{xo + l)|*o) - exp(-i2^Q,^+i)| (102) 

is bounded by J/7 times a function decaying superpolynomially in L plus poly(_L) J/A^ times a function Thus, for 
any xq we have 

Corollary 6. Assume that the system has an {L/8,t) unique bulk state, with t greater than or equal to l/poly(L). 
Assume that the system has a mobility gap with tmax superpolynomially large in L. Then, 



{\\Wiixo)\-^o) - zi\^o)\) (103) 
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is bounded by J/7 times a function decaying superpolynomially in L plus poly(L)J/A^ times a function decaying 
superpolynomially in X^t^ax and also 

{^o\Wiixo)\^o) (104) 

is superpolynomially close to 

xq 

zexp{-i2nJ2Qx), (105) 

x=0 

for some z which is independent of Xq. 

Thus, unless is an integer for all x, we have defined a flux operator which has an expectation value which 
depends on the particular line xq we choose. If we considered the case of a translation invariant system, this would 
be a contradiction, and would prove that there is not a mobility gap and a unique ground state (this is how the 
proof of the higher-dimensional Lieb-Schultz-Mattis theorem goes). In this case, we simply identify that there is a 
position-dependent expectation value of a flux insertion operator. 



LIEB-ROBINSON BOUNDS FOR QUASI-ADIABATIC CONTINUATION 

The previous section relied on the fact that quasi-adiabatic evolution can be approximated by a sum of local 
operators, and hence the evolution D could be approximated to superpolynomnial accuracy by a sum of I?^^-* 
We now want to consider a stronger property. Suppose we have an operator, such as V which can be approximated 
by a sum of local operators. For example, suppose V can be approximated to superpolynomial accuracy, in I, by a 
sum of operators supported on sets Z of diameter at most I. Then, if we had an the additional bound on 
then we would have a Lieb-Robinson bound for the unitary evolution 

Us^S'c^p{i [ ds'V,,}. (106) 
Jq 

In the next two sections, on Hall conductance, we will rely on an assumption of a Lieb-Robinson bound for the 
quasi-adiabatic evolution operator. This bound is fairly immediate to prove if we consider the slightly simpler case of 
a mobility gap rather than a spectral gap. For example, consider the operator defined in Eq. (|3|. For simplicity, we 
can use an exact quasi-adiabatic evolution operator. Then, we have a bound on the operator ^{Hs, dsHz{s)) bounded 
in norm by a constant times {l/^)\\dsHz{s)\\, and the operator also decays superpolyomially in space. Hence, we 
have a Lieb-Robinson bound(see the appendix for more discussion of this case). 

However, in the case of a corrected quasi-adiabatic continuation operator, we also have to worry about the contribu- 
tion from states below the mobility gap. In this case, the contribution of these states, the term J F{\'^t)W!:^°Q{0,t)dt 
in the corrected quasi-adiabatic evolution operator, is local as in lemma ([5|, but the bound on the norm is quite 
weak. In particular, the norm may scale with 1/A^, and hence may scale polynomially with L. This makes it difficult 
(perhaps impossible) to directly prove the Lieb-Robinson bound for corrected quasi-adiabatic evolution directly from 
the assumption of a mobility gap, so we will need one additional assumption. 

In this section, we will define the property of a Lieb-Robinson bound for corrected quasi-adiabatic evolution that we 
need in the next two sections. We also present one simple assumption on the local density of states under which this 
Lieb-Robinson bound can be derived. Then, the results in the next two sections will depend either on the assumption 
of the Lieb-Robinson bound for corrected quasi-adiabatic evolution, or on the (fairly mild) assumption on the density 
of states. 

We will consider parameter-dependent Hamiltonians in the next two sections where flux is inserted along lines. 
We will consider, however, more lines than in the previous section. In the case of a torus, we will have two lines 
describing flux inserted in one direction (the "horizontal" direction of the torus) and two other lines describing flux 
inserted in the other direction (the "vertical" direction of the torus) . The reason we have two lines in each direction 
is similar to the case in the above section: we use the fact that if we insert opposite flux on two different lines then 
the Hamiltonian is only changed by a unitary transformation. The reason we need two different directions is that the 
Hall conductance is equal to the curvature when transported around an infinitesimal loop in fiux space. 

We will need the following Lieb-Robinson bound: 
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Definition 7. Consider a particular parameter dependent Hamiltonian Hs- Consider the quasi- adiabatic evolution 
operator, at s = 0, and corresponding unitary Us- Then, we say that two sets A, B are separated if, for s of 
order unity, we have that for any operator Oa supported on A and Ob supported on B that \\[U sO aUI , O b\\\ and 
\\[UjOAUs,OB]\\ are both bounded by \ A\\\0a\\\\0b\\ times a quantity superpolynomially small in dist(A,i?). 

This Lieb-Robinson bound can be proven, as explained above, for a system with a spectral gap. If there is a 
mobility gap, it can be proven under the following assumption (with slight modification in the definition of the 
corrected quasi-adiabatic continuation operator) . We will consider later parameter-dependent Hamiltonians such that 
dsHs = J2z ^sHz{s) is supported on the sets of sites within distance R of two vertical lines (the solid and dashed 
vertical lines in Figs. 2 and 3). The derivatives dsHz{s) will be non- vanishing only for sets Z which are within distance 
R of one of these lines. We will be interested in sets A and B which are connected by part of one of these lines. 

Consider a given set Z. The operator dsHz{s) for that set has matrix elements between the ground state and 
various excited states. We can define a density of states, pz{E) by 

Pz[E) = {l-PE)(dsHz{s))\^o)\ (107) 



where Pe projects onto states with energy Eq -\- E oi more and the partial derivatives are taken at s = (later when 
we use this density of states, we will always be considering Hamiltonians which are unitarily equivalent for different 
s, so the density of states pz{E) will be independent of s). 

While it is expected that a disordered system will have states with energy of order l/V, we expect that only for 
Z close to certain points will dsHz{s) have non-negligible matrix elements to these states. In contrast, most Z are 
expected to have the property that dsHz{s) will only produce non-negligible matrix elements from the ground state 
to excited states with energies of order unity. In fact, we will require even weaker conditions than that. We will allow 
a typical Z to have non-negligible matrix elements to energy which are of order 1/L", for a < 1. 

Suppose there is an energy A ^ "y/L" such that the following property holds. We define SCOW to be the set of Z 
for which p(A) is bounded by Pmax times ||9s-ff2(s)||, where Pmax is a quantity which is superpolynomially small in 
L. We define J- AST to be the remaining Z. Then, we define the corrected quasi-adiabatic continuation operator as 



I?,=^2?f (108) 



as in Eq. (|2]J. However, for Z G SCOW, we define Vf with the cutoff A< = A, while for Z G TAST wc use the 

cutoff A< = Xjnin- 

Then, this operator continues to approximate the exact evolution up to superpolynomially small error bounded 
by {pmax/ X^)J2z£SCOW \\(^sHzis)\\, with the additional error due to corrections from states below energy A. 

Now, this operator will separate sets A and B under mild assumptions on the density of states pz{E). The 
operators Vf for Z e J- AST will have large norm, but will decay exponentially in space. Assume Xmm > po\y{l/L). 
If we can find a segment of the line of length separating sets A, B as shown in Fig. 1. which scales as /j, for some 
with a < /3 < 1, such that all Z in that segment are in SCOW, then we will have the desired Lieb-Robinson bound: 
the Z e SLOW will give operators 2?f bounded in norm by 0{{L°' /j)\\dsHz{s)\\), and so have a Lieb-Robinson 
velocity, the Z e FAST will have some effect on the dynamics in this segment, due to long-distance tails of dsHz{s) 
(i.e., even is Z g FAST, the operator dsHz{s) has support in SCOW)] however, this produces only corrections of 
order poly(L) times a quantity decaying exponentially in L". Hence, we will have the desired separation. 



HALL CONDUCTANCE ON A TORUS 



We now prove Hall conductance on a torus. For a site i, we define it to have x and y coordinates x{i) and y{i). 
We consider a Hamiltonian with a mobility gap, with t^ax superpolynomially large in L. We now define We define a 
parameter-dependent Hamiltonian H{6i, 02, 4>i, 4>2) as follows. We pick two vertical lines at xi, X2 and two horizontal 
lines 2/1 , y2 to insert flux. Let Qx be defined by 

Xi'Cx{i)'Cx2 

Qx= (109) 

i 

where x{i) is the ai-coordinate of site i and 

yi<y(i)<y-2 

Qy^ ^ q^. (110) 
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(0,0) 



x=L 



FIG. 1: Line illustrating an example of separation on the torus. Sets A and B are marked by circles; there are periodic boundary 
conditions so B is a contiguous set. The flux is inserted along the dashed line connecting them. The short solid intervals mark 
sets SLOW separating A and B as described. 



Definition 8. Let H be any operator which can be written as H = Hz with the Hz supported on a set Z of 
diameter less than L/2. Assume that all the sets Z are squares. Then, each such Hz intersects at most one of the 
lines X — Xi or x — X2 and at most one of the lines y — yi or y — yi. Then, define the twisted operator 



where Hz{0i, O2, 4>i, 4>2) is defined as follows. If the set Z intersects the vertical line x = xi, then Hz{Oi, 02, </>!, 02) = 
exp(i6'i(5x)-ffz(0, 0, 01, (^2) exp(— i^iQx); if the set Z intersects the vertical line x = X2, then Hz{0i,02,4'i,4>2) = 
exp(— «02<9x)-f^z(O, 0, 01, 02) exp(i6'2(9j(:); otherwise Hz{0i,02i4>ii4'2) — ^z(0, 0, 0i, 02)- If the set Z intersects the 
horizontal line y = yi, then Hz(0, 0, 0i, 02) — exp{i(j)iQY)Hz exp(— i0i(5y); if the set Z intersects the horizontal line 
y = 2/2, then ^^(0, 0, 0i, 02) = exp(-i02Qy)i?z exp(i02Qy); otherwise iJz(0, 0, 0i, 02) = Hz- 

Note that we chose the sets Z to be squares so that they would be contiguous sets; thus, if Z was close to a line 
X = xi and contained some points with x < xi and some with x > xi then Z will intersect the line x = xi, and 
similarly for the other three lines. This definition defines our parameter-dependent Hamiltonian, but we will also use 
it later for other operators. 

In this section, we prove the following: 

Ttieorem 2. Consider a Hamiltonian with a mobility gap, with t„iax superpolynomially large in L, and Xmin greater 
than or equal to l/poly(i). Suppose we can find lines Xi,X2,yi,y2 with \xi — X2\ = L/2 and \yi — 2/2I — i/2 such that 
the following holds. Consider any pair of lines. Let A be the set of points within distance L/8 of the intersection of that 
pair of lines. Consider any other pair of lines. Let B be the set of points within distance L/8 of the intersection of that 
pair of lines. Consider quasi- adiabatic evolution under any of the four different quasi- adiabatic evolution operators 
256»,i^2?6i,2 250,1 250,2 figji^gfi below. Suppose that A and B are separated under quasi- adiabatic evolution under all 
four such operators. Then, the Hall conductance is quantized to n(e^//i), for some integer n, up to an error which is 
superpolynomially small in L. 

Note that by assumption the Hamiltonian has an a (i?, r) unique ground state, with r greater than or equal to 
l/poly(L). We only sketch the proof, since it essentially follows [S] (this proof has some similarity with the Chern 
number approach |20| but avoids any averaging assumptions). We pick a parameter r which is superpolynomially small 
in L (the choice of r will be given later), with 27r/r — N for some integer N . Following j^, we define a set of set of 
N"^ different closed paths in parameter space. Each path keeps 62 — <j>2 — throughout. Each path is defined by a 
given pair of numbers {0x,Oy) with 0^ — mr and 0y = nr for a pair of integers m,n in the range 0, ...,A'^ — 1. The 
paths start at 0i = 0i = 0. Then, we move to 0i = 0, 0i = 0y. Then we move to 0i = 0x, 0i = 0y Then we move 
around a small square loop of size r at 0^, 0y as follows: we move to 0i = 0x f, 0i = 0y, to 0i — 0x + 0i = + r, 



We fix \xi - .T2I = \yi - y2| = L/2. 




(Ill) 



z 
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y=L 



(0,0) x=L 

FIG. 2: Lines illustrating how the twists are defined on the torus. The twists 9i, (j>\ affect interactions close to the vertical and 
horizontal solid lines, respectively, while the twists 62 , 02 affect interactions close to the vertical and horizontal dashed lines. 

^1 = ^x, 01 = + T, and back to 9i = 9x,(t>i = Oy Finally, we return to the origin by moving to 9i = 0, 0i = 9y and 
then to 9q — — 0. 

Thus, each path consists of 8 different line segments. We break it into three distinct parts. First, V{m,n) moving 
from 0,0 to to 9i =0,01 = 9y and then to 9i = 9x,(f>i — 9y. Then, Vo(m,n) moving around the small square 
loop. Finally, V^{m,n) returning to the origin. We decompose the first part, V{in,n) — U2{'m,n)Ui{m,n), with 
Ui{m,n),U2{m,n) corresponding to the two segments of that part. We also define unitaries Ui{m,n),1J2{m,n) 
corresponding to evolution along the same path of 0i,0i but with 6*2 = —9i and 02 = —'Pi- We also decompose the 
motion around the square loop as V0{m,n) = Ui){rn,n)UL{m,n)Uijifn,n)Uj^{'m,n) corresponding to moving around 
the four different sides of the square (we use subscripts D, L, U, R indicating the the motion is first to the right, then 
up, then to the left, then down). We also define {/^(m, n), {/^(m, n), Uuim, n), [/jj(m, n) corresponding to evolution 
around the same square but with 02 = —9i and 02 = — 0i. 

Given these paths, we have to pick a quasi-adiabatic evolution operator to generate the various unitaries U and 
U. We define 'D^{9,(t>) to be the quasi-adiabatic continuation operator which produces an infinitesimal change in 9 
starting at 9i — —02 = 9 and 0i — —02 — 0. That is, for Hg = H{9 + s, —9 — s, 0, —0), wc define T>^{9, 0) = "Dg at 
s=0. We define 'D'^{9, 0) to produce an infinitesimal change in 0, again starting at 9i = —02 = 9 and 0i = —02 = 0. 
We use these operators 'D^{9,(j)),'D'^{9,(j)) to generate the unitaries Ui{m,n),U2{rn,n),UD,L,u,B.{'rn,n). 

Let us consider the operator 1)^(0,0). For 9i = —02 = s and 0i = —02 — 0, we have that dsH{s) — J2z ^sHz{s), 
and dgHz{s) is nonvanishing if Z is within distance R of the line at x = xi or if Z is within distance R of the line 
at x = X2. Let 0'^^\s) denote the sum of terms in dsH{s) near the line at x = xi and let O'-^^(s) denote the sum 
of terms in dsH{s) near the line at x = X2, so that Vg = 'D{Hs,0'-^^s)) +'D{Hs,0^'^\s)). Under the assumption of 
a mobility gap, because of the superpolynomial localizability property, up to superpolynomially small error, we can 
approximate the operators 2?(-ffs, O^^-*) and 'D{Hs, O'^^) by operators I?*'^^) and I?*'^^) supported within distance less 
than L/8 of the respective lines x = and x = L/2 and which are a sum of squares of operators on squares of size at 
most L/2. Then, we can apply the definition Q to define 

I?«'i(0,0)=I?^^«(0,O,0,O), (112) 

and 

2?'''2(0, 0) = 2?'''(2)(o, -0, 0, 0). (113) 

Note that 

I?'^(0,O) - (2?^'i(0,O) +2?'''2(0^O)) (114) 

is superpolynomially small. We define V^'^^^ and V^'^^^f in an analogous way by decomposing the operator I?"^(0,0), 
and then we have that 



2?^(0, 0) - (2?^'i(0, 0) + V^'\0, 0)) 



(115) 
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is superpolynomially small. We will use these operators V^'^{6,(j)) and V^'^{9,(j)) to generate the unitaries 
C/i(m, n),U2{m, n), UD,L,u,R{m, n). 

The reader will notice one fact about this choice of quasi-adiabatic evolution operators: they are all defined in terms 
of operators at = (/> = by applying a twist. This implies that the separation assumption can be defined solely in 
terms of the operators T) d.t 6,<j> — 0. For example, 

S'exp{iJ^ ds'V'^'\s,0)} = exp(iQxs')S' exp{i ds'(-Qx + 2?^'^(0, 0))} exp(-iQxs')- (116) 

The proof rests on the following four facts. We sketch each in turn. First, for r sufficiently small, the evolution 
around the path U (0, 0) returns the ground state to the ground state up to superpolynomial small error and up to 
some phase which is that is, up to superpolynomially small error, equal to r"^ /2it times the Hall conductance in units 
of e^/h, plus corrections of order r^. That is, we claim that 

Lemma 11. The quantity 

|(1'o,f/(0,0)*o)| (117) 
is superpolynomially close to unity for any r of order unity. Also, the quantity 

|(vI/o,C/(0,0)*o)-exp(ir2^4)l (118) 
is bounded by terms of order poly {L)0{r^) plus terms which are superpolynomially close to unity. 



Proof. Eq. (1181 follows from a power series expansion of the quasi-adiabatic continuation operator. The terms of 
order in the expectation value in Eq. (118) arise from approximating, along each leg of the square, the evolution 



U, = S' exp{i [ Vsds'} = 1 + isVo 
Jo 



+ 0(s2). (119) 



The quasi-adiabatic evolution operator i'Do acting on the ground state is superpolynomially close to the derivative 
of the ground state with respect to the parameter s. Thus, to order r^, the expectation value (^E'o, f/(0, 0)^*0) is 
superpolynomially close to 1 + r'^2lm{{d^'i> {6,(1)), dg'if {9, ((>))), where '^{6,(j)) is the ground state of H{9,0,(t),Q), and 
partial derivatives are taken a.t 9 — (f> = 0. However, this quantity 2lm{{d^'^{9,(l)),dg'i'{9,4>))) is related to the Hall 
conductance by the Kubo formula, so Eq. (118) follows to order r^. Using a bound on the norm of V, we bound the 
terms of order by poly(L). 

To show Eq. (117), we show that the state U{0, 0)'^o has an energy which is superpolynomially close to Eq. Then, 
using the assumed \mim Eq- (117) will follow. We will estimate the energy by estimating the expectation value of 
Hz for each Z in the state [/(0,0)\I'o- 

Suppose first that Z is far from the point x = Xi,y = yi where two lines intersect (we use the word "far" to indicate 
something is a distance of order CL from something, and near otherwise). Note that [/(0,0) = Vo(0,0). We have 
Vo(0,0) = C/i5(0,0)C/l(0,0)C/,7(0,0)C//?(0,0). For simplicity, in this lemma we write Ul = Ul{0,0),Ur = Ur{0,0), 
and so on. We claim that the operator V(3(0, 0) can be approximated to superpolynomial accuracy by an operator 
supported near x = xi,y = yi using the separation assumption and superpolynomial localizability property. So, 
if Z is far, from x = xi,y — yi, then || [i/^, t/(0, 0)] || is superpolynomially small, so {'i/o,U{0,Oy HzU{0,0)'^o) is 
superpolynomially close to (^'o, i?z^'o)- 

We now show that V0(O,O) can indeed be approximated by an operator supported near x — xi,y — yi. We will 
make repeated use of the basic identity for any two Hermitian operators a^, b^ that, defining 



rs2 

s^^s2 = 5'exp{i / ds'6, 

J si 



} (120) 



we have 



5'exp{i / ds'os'} = S' exp{i / ds'us' ,s{as' ~ bs')ul, ^}uo,s (121) 
JO Jo 

= Mo,s5'exp{z / ds'ul, Q{as' ~ bs')us'fi}- 

JO 
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Decomposing Vo(0, 0), we have 

UdUlUuUr = UoUuS'eM-iUl f ds'V^'^r - s',r)Uu}UR (122) 

Jo 

So, it suffices to show that UdUu can be approximated by an operator supported near x = xi,y = yi and similarly 
for S'e^p{-iUl ds'V'^^^r - s,r)Uu}UR. We have 

nr nr 

UdUu = S' exp{-i ds'V'>''^{0,r~s')}S'exp{i ds'V^^^ir, s')}. (123) 

Define ^(s') = V'''-\0,s') - V'l'-^{r,s'). Define Ws^^s^ = 5'exp{i /^'^ ds'X>"*'i(r, s')}, so that Uu = M^o,r. Then, 

UdUu = S' exp{-^ / ds'W,,^rA[r - s')wl ,,}. (124) 
Jo 

Note that A{s') is can be approximated by an operator supported near x = xi,y = yi using the supcrpolynomially 
localizability property, and so using the separation assumption, Ws\rA{r — s')Wl, ^ can also be approximated by such 
an operator to superpolynomial accuracy. Thus, UdUjj can be approximated by such an operator. We next consider 
S' exp{-iC/^; /g ds'V'^-^r - s', r)Uu)UR. Define 

B{s) = uljV'^'^r - s,r)Uu -V^'^ir - s,r), (125) 
and define Y,^,^^ = S' exp{-i J^^ ds'V'^'^s - s',r)} so that Ul = Yo^r- Then, 



S'exp{-iUl / ds'V'''\r-s,r)Uu}UR = S'exp{-i / ds'Ys^ ,rB{s')Y^, ^}UlUr. 
Jq Jo 



(126) 



Note that B(s') is can be approximated by an operator supported near x = xi,y = yi using the supcrpolynomially 
localizability property, and so using the separation assumption, Yg'.,._B(s')yJ, ^ can also be approximated by such an 

operator to superpolynomial accuracy. Thus, iS'exp{— i ds'Yg/ ,rB{s')Y^, ^} can also be approximated by such an 
operator. Finally, one may show that UlUr can be approximated by such an operator in a similar way as to the proof 
that UdUu could be approximated by such an operator. 

Now, suppose that Z is close to both y — yi and x — Xi. Consider the expectation value 
i^OjUj^UljUlujjHzUDULUuUR'^o) ■ Note that uj^HzUo is superpolynomially close to Uj^HzUo- So, we may 
consider {'^q^U^uIiU'IU^iHzUdUlUuUr'^q). Using the separation assumption, uj^HzUo has small commu- 
tator with 2?^'^. So, u\u]-,HzUdUl is superpolynomially close to U\u]-,HzUdUl (this calculation uses the 
identity ( 121[ ) and similar decomposition of unitaries into ordered exponentials as the above case, so we omit 



the details). Repeating this, we find that {'^q.uI^uIjUIu^iHzUdUlUuUr'^q) is superpolynomially close to 
{^a,UlfUljU\ij]^HzUDULUuUR^{^,). However, UdUlUvUr^o) is superpolynomially close to 'So times a phase. 
Therefore, {'^q,u\iUIjU\u\)HzUdUlUuUr'^q) is superpolynomially close to {^a.Hz'^o) ■ This completes the possi- 
ble cases. □ 



Let us discuss one point in the above proof. This approach used to show Eq. (117) is different than the approach 
used in [8], where it was simply noted that quasi-adiabatic evolution around a closed loop in a gapped region of 
parameter space approximates adiabatic evolution and hence returns the ground state to the ground state up to small 
error; we use this different approach because here we do not want to rely on any assumption of a mobility gap for 
parameter values other than 9 — <f> = Q. In [S], since r was taken superpolynomially small, evolution within the region 
of parameter space with \9x\, \0y \ < r preserved the spectral gap assumed in that paper. In this paper, we could try 
to follow that idea and show that a mobility gap at = (p = implies a mobility gap for superpolynomially small 
change in the Hamiltonian; this may be possible, but since an alternate method is available we preferred not to do 
this. 

The second property is that the evolution around all loops returns the ground state to the ground state and produces 
the same expectation value up to superpolynomially small error. That is. 

Lemma 12. Define 



z{m,n) = {•^o,U{m,n)^o). 



(127) 



23 



Then, 



is superpolynomially small. 



zofi\ (128) 



Proof. To prove this, note that the expectation vahie is equal to {'^Q\V{m^n)^V(^{m^n)V{m,n)\'^i^) . The op- 
erator can be approximated to superpolynomial accuracy by an operator supported near x = Xi,y — 
yi. So, U2{m,nYVQ{m,n)U2{m,n) is superpolynomially close to U2{m,nyVQ{rn,n)U2{rn,n). Again us- 
ing the Lieb-Robinson bounds, Ui{m,nyU2{m,nyV(^{m,n)U2('m,n)Ui{rn,n) is superpolynomially close to 
Ui{m,n)W2{m,nyVQ{m,n)U2{m,n)Ui{rn,n). However, U2{m,n)Ui{m,n)\'^Q) = exp{iQxdi) exp{iQY02)\'^o) 
and 1^0(771, n) is superpolynomially close to ex.p{iQxOi) exp{iQYd2)V0{m,n) exp{—iQx0i) sx.p{iQY02)- So, 
{^o\U{m,n)\^/Q) is superpolynomially close to (5'o|J7(0, 0)|*o)- □ 

The third property is that the product (u{N ~ 1, N - 1)U{N -2, N ~ 1)...U{0, N - 1)J (u{N - 1, N -2)....U{0, N - 

2)j ...(u{N — l,O...C/(0,0)^ is exactly equal to the unitary operator Ubig corresponding to quasi-adiabatic evolution 
around a big loop of size 2tt: moving from 6'i = 0i = to 6*1 = 27t, (pi = 0, to 0i — (pi = 27r to 9i — 0,(pi = 2tt, to 
9i ^ (pi ~ 0. This is the evolution decomposition in [8j. 

The fourth property is that the phase due to evolution around the big loop is equal to unity up to superpolynomially 
small error. To prove this, we need to show that evolution along any of the four sides of the big loop (for example, 
from 6'i = 01 = to ^1 = 2n, (pi = 0) send the ground state to the ground state up to a phase, up to superpolynomially 
small error. Then, the phase will cancel between opposite sides of the loop. However, the fact that evolution along a 
single side sends the ground state to the ground state up to a phase and up to superpolynomially small error follows 
from the same argument as in the Lieb-Schultz-Mattis proof previously. Thus, 

{^o,UMg'^o) (129) 

is superpolynomially close to unity. 

We then approximate the phase of evolution around the large loop by the product of phases for evolution around 
each small loop, following the Appendix in [8]. That is, we approximate: 

The error in the second approximation is due to the fact that each small loop has some "leakage" outside the 
ground state. That is, U{m,n)^Q has some component perpendicular to the ground state. Define E^.n to be the 
1(1 - Po)U{m,n)'^o\^ = 1 - \zm,n\^- Then, the difference 

Y^V^n. (131) 



However, by lemma (11), the quantity £'0,0 is superpolynomially small and so by lemma (12), the quantity E^^n is 
superpolynomially small for all m, n. The number of small loops scales as l/r^. Hence, we can pick an r which is 



superpolynomially large such that the difference in Eq. ( 131 ) is negligible. Thus, |1 — Hm n ^rn.nl is superpolynoomially 
small. Then, since the phase is almost the same for every small loops, the phase around the loop with m = n = is 
close to 2TTn/N — 27r(r/27r)^ for some integer n. Hence the Hall conductance is close to ne^ /h. 

Note that we used three properties of these operators 'D^{6,(p) and 'D'^{6,(p). First, that each operator can be 
approximated, to superpolynomial accuracy, by a sum of operator operators supported close to the respective lines 
X — xi,X2 or y — 2/1,2/2, each of which is approximated by a sum of operators on squares of increasing size but 
superpolynomially decaying strength. Second, the separation assumption. Third, evolution of the ground state under 
T)^ or T)'^ approximates, to superpolynomial accuracy, the adiabatic evolution of the ground state, and hence gives a 
state which is is superpolynomially close to exp(?(5j!c0) exp(i(5y(/))\E'o. Thus, any operators which satisfy these three 
properties would suffice. 



24 



HALL CONDUCTANCE ON AN ANNULUS 



We now consider open boundary conditions. Specifically, we consider a system on an annulus. In this case, a system 
will not have a gap: there will be gaplcss edge modes. However, there is still a bulk gap. We need to define this. 
We let T and O denote the sites on the inner and outer edges of the annulus. We let B denote the sites .s such that 
dist(s,/) > L/3, dist(s, O) > I//3. That is, B is a ring around the inside of the annulus, which we call the "bulk". 
The constant factor 1/3 is not particularly important; we simply want a constant less than 1/2 (so that the widht of 
the ring scales linearly in L) and greater than (so that the distance from the edges also scales linearly with L). We 
use X and y coordinates to label sites, and let the line ?/ = correspond to halfway between the inner and outer rings 
of the annulus, so that B includes points with — i/6 < y < L/6. 

Definition 9. We say that a system has a bulk mobility gap 7 and localization length ^ a,nd localization constant 
cioc up to time t„iax if, for any operator O supported on set A, there exists an operator Wl°Q{0,t) with the following 
properties. First, for any I, W^°g(0, t) can he approximated by a sum of operators Pi + Po + Pbuik, wth Pj is supported 
on bi{I), Po is supported on bi{0), and Pbuik is supported on bi{A), up to an error in operator norm bounded by 

Q„eexp(-//C)||W^^,G(0)|| +max|^l>^|G(w)|||0|| (132) 

with 

\\Pi\\ < ||O||Q„eexp(-dist(A^)/0. (133) 



IIPoll < ||0||Qocexp(-dist(A,0)/e). (134) 

Second, we require that the state produced by acting with Wj^G{0){t) = exp(iffoi)M^7,G(0) exp(— iifoi) on the ground 
state is equal to the state produced by acting with an operator W}^"Q{0,t) on the ground state, i.e., 

VF.,,g(O)(f)*0 = W^:G(0,i)*o. (135) 

Third, we have 

\\wi;;G{o,m<\\w,,Gm. (136) 

Definition 10. We say that a Hamiltonian H has an {I, r) unique bulk ground state if the following holds. Given 
any density matrix p such that, for all sets A c bi{B) with diam(A) < I we have 

\\Tr^{p-Po)\\i<e, (137) 

then 

\\TT^{p-Po)\\i<rV~e (138) 

The definition of a unique bulk state is certainly necessary. Let us explain why. Consider a fractional Hall system 
defined on an annulus. The system has a bulk gap. However, the system does not have a unique bulk state according 
to our definition: physically, one expects q different bulk states in a l/q Laughlin wavefunction. Such a system will 
also not have an integer Hall conductance. Thus, to prove integer quantization we need to make the assumption of a 
unique bulk state. 

We now define twisted boundary conditions. We pick two vertical lines at Xi,X2 with \xi — X2\ = L/2 and a single 
horizontal line at y = 0. Let Qx be defined by 

Xl<x{i)<X2 

Qx= (139) 



where x{i) is the f-coordinate of site i and 



yi>o 

Qy=Y1 1i- (140) 




FIG. 3: Lines illustrating how the twists are defined on the annulus. The twists 61,62 affect interactions close to the vertical 
solid and dashed lines, respectively, while the twist affects interactions close to the horizontal dashed line. Regions X and O 
are the top and bottom rectangles filled with angled lines, respectively. 



Definition 11. Let H be any operator which can be written as H = Hz with the Hz supported on a set Z of 
diameter less than L/2. Assume that all the sets Z are squares. Then, each such Hz intersects at most one of the 
lines X = xi or x = X2 and at most one of the lines y = yi or y = y\. Then, define the twisted operator 

H{eu02,(l>) (141) 

as follows. Let 

if(^i,^2,<^) = $^ifz(^i,e2,0), (142) 

z 

where ffz(6'i, 6*2, </>) is defined as follows. If the set Z intersects the vertical line x = Xi, then Hz{0i,92,(p) = 
exp{i9iQx)Hz{0,0,(j))exp{—i9iQx); if the set Z intersects the vertical line x = X2, then Hz{0i,62,4'i) = 
exp{—i92Qx)Hz{0,0,(j))exTp{i62Qx); otherwise Hz{0i,62,(f)i,^2) = Hz{0,0,^). Finally, define Hz{0,0,(t)) = 
exp{i(pQY)Hz exp(-i^Qy). 

Note that there is only one horizontal line for the annulus, while we had two for the torus. 

In this section, we prove the following: 

Theorem 3. Consider a Hamiltonian with a bulk mobility gap, with tmax superpolynomially large in L, and Xmin 
greater than or equal to l/poly(Z/), and with a {cqL,t) unique bulk state, with r greater than or equal to l/poly(i), 
for sufficiently small Cq. Suppose we can find lines Xi, X2 with \xi — X2I = L/2 such that the following holds. Let A be 
the set of points within distance L/6 of the intersection of the line x = X2 with y = 0. Let B be the union of the set 
of points within distance L/6 of the other line and the set of points within distance L/6 of either edge of the annulus. 
Consider quasi- adiabatic evolution under the quasi- adiabatic evolution operators defined below. Suppose that A and 
B are separated under quasi- adiabatic evolution under all such operators. Then, the Hall conductance is quantized to 
n{e^ /h), for some integer n, up to an error which is superpolynomially small in L. 

Note that we allow an operator on the line x = x\to evolve under quasi-adiabatic evolution to an operator which 
cannot be approximated by an operator supported near x = Xi. We allow the operator to spread out along the line, 
and along the boundary of the sample, having a large effect on sites along the boundary, but we require that it not 
penetrate into the bulk near the line x = X2. This is needed for the second property below. 

The proof essentially follows the above proof. We again define a set of set of N"^ different closed paths in parameter 
space. Each path keeps 92 = (j) = throughout. Each path is defined by a given pair of numbers {O^, 9y) with 9x = mr 
and 6y = nr for a pair of integers m, n in the range 0, A'' — 1. The paths start at 9i = (p = 0. Then, we move to 
9i = 0,(j) = 9y. Then we move to — 9^, ((> = 9y. Then we move around a small square loop of size r at 9^, 9y as 
follows: we move to ^1 = 6x + r,4> = 9y, to 9i = 9x+r,<f) = 9y-\-r, to 9\ = 9x,(p = 9y + r, and back to 9\ = 9x,4' = 9y. 
Finally, we return to the origin by moving to = 0, </> = Oy and then to 9q = (f) = Q. 
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Thus, each path consists of 8 different line segments. We break it into three distinct parts. First, V{m,n) moving 
from 0,0 to to 0i — 0,(j) = Oy. Then we move to Oi = 9xt4> = dy Then, Vo(m,n) moving around the small 
square loop. Finally, V^{m,n) returning to the origin. We write the first part, V{m,n) = U2{m,n)Ui{'m,n), with 
Ui{m,n),U2{'rn,n) corresponding to the two segments of that part. We also define unitary Ui{m^n) and U2{m^n) 
corresponding to evolution along the same path of 0i, but with 62 — —01- 

Given these paths, we have to pick a quasi-adiabatic evolution operator to generate the various unitaries U and 
U. We define 0) to be the quasi-adiabatic continuation operator which produces an infinitesimal change in 9 

starting at 9i — —92 = 9 and the given </). That is, for Hg = H{9 + s,—9 — s,(f>), we define 'D^{9,(j)) — Vg at s=0. 
We define 'D'^{9, (p) to produce an infinitesimal change in </>, again starting at 9i — — 6*2 = 9 and the given (f). We use 
these operators {9 , 4>) {9 , <j)) to generate the unitaries IJi{rn,n),U2{m,n),Uo,L,u,F!.{^Tn). Let us consider the 
operator I?^(0,0). For 9i = —92 = s and = 0, we have that dsH{s) = dsHz{s), and dsHz{s) is nonvanishing 
if Z is within distance R of the line aX x = xi or if Z is within distance R of the line at a; = X2. Let O'^-'^^(s) denote 
the sum of terms in dsH{s) near the line at x — xi and let 0'^'^\s) denote the sum of terms in dsH{s) near the 
line at a; = X2, so that Vg = 'D{Hs,0^^\s)) +'D{Hs,0'^-^\s)). Under the assumption of a mobility gap, because of 
the superpolynomial localizability property, up to superpolynomially small error, we can approximate the operators 
'D{Hs,0^^'') and by operators V^P and v'^P supported within distance less than L/8 of the respective 

lines x — Q and a; = i/2 and which are a sum of squares of operators on squares of size at most L/2. Then, we can 
apply the definition Q to define 

2?'''i(0,0)=I?W (0,0,0), (143) 

and 

V'^\9,cj,)=Vf\Q,-9,(ly). (144) 

Note that 

V%9, 0) - {v'^^^{9, 0) + V^'^{9, 0)) (145) 

is superpolynomially small. Similarly, we define T)'^'^ as follows: we approximate I?'^(0,0) by an operator calDg which 
is a sum of operators supported on squares of size at most L/2. Then, we apply definition ([s]) to define 

V^^\9,,p)^V,{9,{),<jy). (146) 



Note that I?*(6', 0) is superpolynomially close to 'Ds{9,—9,(f>), so at 7^ 0, the operators I?"^'^(6', 0) and 
are not necessarily superpolynomially close to each other. We will use these operators 'D^'^{9,(f)) and 'D^'^{9,(f>) to 
generate the unitaries Ui{m,n)^U2{Tri^n),UD,L.u.R{^in). 

The proof rests on the following four facts. We sketch each in turn. 

First, for r sufficiently small, the evolution around the path U{Q, 0) returns the ground state to the ground state up 
to superpolynomial small error and up to some phase which is that is, up to superpolynomially small error, equal to 



r /27r times the Hall conductance in units of e //i, plus corrections of order r . That is, similar to lemma (11) above, 
we claim that 

Lemma 13. The quantity 

|(^'o,C/(0,0)*o)| (147) 
is superpolynomially close to unity for any r. Also, the quantity 

|(vI/o,C/(0,0)*o)-exp(zr2^A)| (148) 

is bounded by terms of order poly {L)0{r^) plus terms which are superpolynomially close to unity. 
The proof of this is analogous to the torus case. 

The second property is that the evolution around all loops returns the ground state to the ground state and 
produces the same phase up to superpolynomially small error. To prove this, note that the desired expecta- 
tion value is {'i>o\V{m,n)^V(^(m,n)V{m,n)\'^o)- The operator V(j can be approximated to superpolynomial ac- 
curacy by an operator supported near x = xi,y = j/i, using the Lieb- Robinson bound for quasi-adiabatic con- 
tinuation. So, U2{m,n)V(^{m,n)U2{m,n) is superpolynomially close to U2{m,nyV(^{m,n)U2(rn,n). However, 
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C/2(to, n)'''/7i(m, n)|^'o) = exp(iQx^i) exp(z(5y(^)|^'o) for 9i = mr,(j) = nr, and l^(m,n) is superpolynomially close 
to exp(«(5x^i) exp(i(5y(^)Vo('Ti, n) exp(— iQjsf 0i) exp(— iQyc/)). So, (^'o|J7(m, n)|^'o) is superpolynomially close to 
(^'o|?7(0, 0)|^'o)- This argument is almost identical to the torus case, except different sets were required to be sepa- 
rated. 

The third property is that the product (u{N -I, N -1)U{N -2, N -1)...U{0, N -1)^ (u{N -1, N ~2)....U{0, N - 

2)j ...(^U (N — 1, O...C/(0, 0)^ is exactly equal to the unitary operator corresponding to quasi-adiabatic evolution around 
a big loop of size 2tt: moving from 6'i = 0i = to 6*1 = 27r, (pi — 0, to 9i — — 27r to 9i = 0,(j)i = 2tt, to 9i = (pi — 0. 
This is the evolution decomposition in [5]. 

The fourth property is that the phase due to evolution around the big loop is equal to unity up to superpolynomially 
small error. To show this, note that evolution along the sides of the loop where (p increases from to 27r, or 
decreases from 27r to zero, is given by an operator that is, up to superpolynomially small error, supported in B. 
Define this operator to be 1/2^^. Acting on the ground state, this operator [/|^^ only produces a phase, since it 
approximates the exact evolution under a unitary transformation of the Hamiltonian. This phase is non-trivial, being 
Z0 = exp(-z27r(3y), where Qy = (*o|Qy|*o)- 

Let f7|jr define the evolution along the side of the big loop where 9 increases from to 2tt. We are interested in 
the action of f/^^ on the state J7|^|5'o). We claim that J7|^^t/|jr|^o) ~ ^0C^27rl*o)j and we will show this in the next 
paragraph. Given that this is true, then the desired fourth property follows: 

(*o|(f/t)^(c/t)Vtc/|j*o) « ^0^0(*o|(t^L)^f^tl*o> = 1. (149) 

To show the claim, it suffices to show that the the reduced density matrix of J/ajrl^o) on B is close to the reduced 
density matrix of I^Pq) on B since ?7|^ is approximated by an operator supported on B. To show this, we use the 
assumption of a unique bulk state: if we can show that the expectation value in the state [/27rl^o) of all operators O 
supported on B supported on sets of diameter at most cqL is close to that in the ground, then we are done. If O is 
distance of order L from x — Xi, then this follows from the locality of the quasi-adiabatic continuation operator. If 
O is close to X — xi, then for sufficientUy small cq, the operator O is supported a distance of order L from x = a;2 by 
assumption. Define 

i?2x = S'exp{i / ds'viP +V^^^} (150) 
Jo 

= C/|,5'exp{z£^d.'(C/f,)^I?(^)C/5}, 

where C/f = 5'exp{z J^^ ds'vl,]^}. Using the assumption of separation, the operator 5'exp|i J^^ ds' (U^,)''V^'^^U^,^ 
can be approximated to superpolynomial accuracy by an operator supported away from the intersection of B with 
X = xi. So, 

{^o\{uLy OUll^o) « (vI'o|i?LOi?2.|^'o). (151) 

However, i?2ir describes quasi-adiabatic continuation along a unitarily equivalent path of Hamiltonians, and i?27r|^o) = 
|4'o). So, this completes the proof that the expectation value of O is almost the same for all operators O of diameter 
at most CqL supported in B. 

Given these four properties, we can complete in the same manner as in the torus case: we approximate the phase of 
evolution around the large loop by the product of phases for evolution around each small loop. The number of small 
loops scales as 1/r^, but we can find an r which is superpolynomially small but for which this leakage is negligible. 
Hence, the product of the phases is superpolynomially close to unity. Then, since the phase is almost the same for 
every small loops, the phase around the loop with m = n = is close to 2Trn/N = 27r(r/27r)^ for some integer n. 
Hence the Hall conductance is close to ne'^/h. 

DISCUSSION 

This paper mostly consisted of a series of definitions. Once the correct definitions are found, the corrected quasi- 
adiabatic continuation operator could be constructed straightforwardly. This allowed us to carry over many of the 
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previous results obtained with these operators in the case of a spectral gap in a straightforward way to the case of a 
mobiUty gap. It is likely that there are other applications of these ideas. 

In the appendix, we presented quasi-adiabatic continuation operators that have very rapid decay in space, simpli- 
fying and tightening previous estimates. 

Acknowledgments I thank C. Nayak and T. J. Osborne for useful discussions on disordered systems. I thank S. 
Bravyi and S. Michalakis for useful discussions on filter functions. 



Optimized Quasi- Adiabatic Continuation 

In this section, we construct an infinitely differentiable function F(u)) such that F{uj) — —1/uj for |a;| > 1 and such 
that the Fourier transform F{t) decays exponentially in a polynomial of t. Our general strategy to construct this 
function is to construct an infinitely differentiable function equal to 1 for \io\ > 1 and vanishing at uj ~ 0, and then 
multiply this function by l/w. 

Then, in the next section we show the application of this function to construct an exact (in that it can be used to 
exactly describe evolution of the ground state of a gapped Hamiltonian under a parameter change) quasi-adiabatic 
operator with good decay properties in space (we present this in detail only for a system with a spectral gap), and 
show Lieb-Robinson bounds for this operator and discuss the best possible bounds. Then, we present some discussion 
of the problem of optimizing over functions. Finally, we present applications in the last section. 

Let f{uj) be even and have the property that /(w) = for = and /(w) = 1 for |a;| > 1. Let f{t) be the Fourier 
transform of /(w). Then, define F{t) by 

Fit) ^'-J duf{u)sigiiit - u), (152) 

where sign(i — u) is the sign function: sign(i — u) — I for t > u, sign(t — u) — — 1 for t < u, and sign(O) — 0. We now 
show the time decay of F{t) and we show that the Fourier transform F{ll!) is equal to — l/w for j^l > 1, as desired. 



Lemma 14. Let F{t) be as defined in 152 Let F{(jj) be the Fourier transform of Fit). Then, 



/•oo 

\F{t)\ < I / f{u)dul (153) 
J\t\ 

and 

F{iu) = — />). (154) 
Proof. Assume, without loss of generality, that t > 0. Then, we have \F{t)\ < \ f{u)du\/2-\- \ f{u)du\/2. Since 



/(O) - 0, we have | f{u)du\ = | /^^ f{u)du\. Thus, \F{t)\ < \ /,°° f{u)d 



We have 



F{u;) ^ ^ J dtexp{iujt) j duf{u)sign{t — u). (155) 

Integrating by parts in t, we have 

F{uj) = J dtexp{iujt) J duf{u)S{t~u) (156) 



-1 
U! 



Note that limf_j.-|-oo ^ J du/(u)sign(t — u)^ = 0, so the contributions to the integration by parts from the upper and 
lower limits of integration vanish. □ 

Thus, given any function f{t) which decays as \ f{t)\ < Cexp(— t"), for some C,a > 0, we can find a filter function 
F{t) which decays as Cexp(— t"), for some other C. Note that given any even function g{t) which has the Fourier 
transform with the property that g(w) = 1 for lu — 0, f{uj) ~ for |a;| > 1, we can define 



f{t)^S{t)~g{t), 



(157) 
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where 6{t) is the Dirac 5-function (since we convolve f{u) against sign(t — u), the resulting F{t) is a function, rather 
than a distribution). 

In the classic paper [l9]. it is shown how to construct such functions g{t) such that 

\git)\<0{cxpHt\ei\m, (158) 
for any monotonically decreasing positive function e{y) such that 

^d, (159) 

1 y 

is convergent. For example, such a function e{y) may be chosen to be 

e(y) = l/log(2 + y)2. (160) 

Thus, our function g{t) has so-called "subexponential decay" [23]. A function f{t) is defined to have subcxponcntial 
decay if, for any a < 1, \f{t)\ < Ca cxp(— i"), for some Cq which depends on a. Thus, 

Corollary 7. There exist odd functions F{t) with subexponential decay and with F(uj) = —1/uj for > 1. In fact, 
the resulting function F{t) from the given e{y) obeys 

\F{t)\ < 0(log(2 + \t\y cxp(-|t|/log(2 + \t\ f)). (161) 



Below Eq. (185), we list even faster possible asymptotic decay. 



Decay of Quasi- Adiabatic Continuation Operator 

Given the rapid decay of F(t) with t, we have a Lieb-Robinson bound for quasi-adiabatic continuation using this 
operator, as we now show. Assume a system has a spectral gap 7, and that Hs = Hs{Z), with Hs{Z) supported 
on sets Z of diameter at most R. Define a quasi-adiabatic evolution operator by 

2?, =^2?f, (162) 
z 

where 

iVf ^ J dtF{jt)exp{iHst)(dsHs{Z)^exp{-iHst). (163) 

If a Hamiltonian Hg has a spectral gap 7 and a ground state 4'o(s), then by construction 

ds^ois) = zI?.*o(s). (164) 

We now consider the decay properties of . Let O be an operator supported on a set which is distance at least I 
from Z. Then, using the Lieb-Robinson bound and a triangle inequality 

||[0,pf]|| < J dt\F{jt)\\\[0,eM^H,t)(^dMZ)^expi~iHM\ (165) 
<([ dt\F{jt)\g{l)+2 [ dt\Fi^t)\) X \\dMZ)\\\\0\\ 

^J\t\<l/VLR J\t\>l/VLR ' 

< E{l)\\dsH,{Z)\\\\Ol 



where we define 



E{1) ^ - ( / du\F{tMl) + 2 / |F(^)|) , (166) 



where we used the change of variables, u — ^t. We define, for use later, F,(fS) = 1. 
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Define bi{Z) to be the set of sites within distance I of Z. Define an approximation to Vf which is supported on 

biiZ) by 



Df ' = j dUU^VfU, (167) 



where the integral ranges over unitaries U which are suspported on the complement of bi{Z) with the Haar measure. 
Using the commutator estimate above, we can estimate the diff'erence WWVfU — || < 'Df]!!- So 

\\Vf'' -Vf\\< E{1). (168) 

Define a sequence of operators Ds{Z,j) by 

D,{Z,0)=Vf'°, (169) 

for j = and by 

D,{Z,j) = Vf'^ -V^^^-\ (170) 
for j > 0. We have the bound for j > 0, that ||£)s(Z, j)|| < 2E{1). Thus, since 

V, = J2Ds{Z,j), (171) 
z,j 

wc have decomposed Vs as a sum of operators on sets of increasing size and decreasing norm as follows. 

For definiteness, let us assume that the initial Hamiltonian H has the property that the sets Z are all given by balls 
of radius R/2 about sets i (let R be even for simplicity). So, we write Hz = Hi, where Z = br/2{{i}), and {i} is the 
set containing just site i. Then, 



Lemma 15. For the given function F{t), 



ieA j>iJ/2 



where the sum is over sites i in the lattice A, and where Ds{i,j) is an operator supported on bj{{i}) and 

\\D,{i,j)\\ < 2E{1 - R/2)\\dsHz{s)\\. (173) 
If a set y = bj{{i}), then we define Ds{Y) = Ds{i,j). It is also worth bounding this as: 
Lemma 16. For any pair of sites i,j, 

^ p,(r)|| <i^(dist(z,j)), (174) 

Y3i,j 

where 

mGA fe>max(dist(i,TO),dist(j,m)) 

where we define 

J' = ma^z{\\dsHz{s)\\). (176) 

We are now ready to derive the Lieb- Robinson bounds for evolution under Vg- Define a unitary Us corresponding 
to evolution with Dg by 

Us=S'e^p{i [ ds'Vs'}. (177) 
Jo 

Let Oa be any operator supported on set A, and let Ob be an operator supported on set B, where B is the set of 
sites which are distance at least I from set A. We wish to bound 

\\[UsOaUIOb]\\. (178) 
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If a set y = bj{{i}), we will write Ds{Y) = Ds{i,j) to save notation in what follows. We use the series expansion 
first derived in pi: 

\\[UsOaUIOb]\\ < 2{2\s\) 11^^(^)11 (179) 

Y-.YnA^HiXnB^ti 

E \mY,)\\\\Ds{Y2)\\ 

Yi,Y2:YnA^(l),YinY2^D,Y2nB^ll) 

+2^^ E \mY^)\\\mY,)\\\mY,)\\ 

Yi,Y2,Y3:YinA^<l),YinY2^il),Y2nY3^il),Y3nB^<ll 

+... 

We may bound this as follows: 

\\[U.OaUIOb]\\ < e BG{i,j), (180) 

where G{i,j) is defined by 

2{2\s\) J2 11^^(^)11 (181) 
Y:ieY,jeY 

,A2\s\r 



^ ||i?.(ri)||p,(r2)|| 



2! 

Yi,Y2:ie>',yinY2#0,iei2 

+2^^ E |!D.(yi)|iiii?.(y2)llP.(i3)ll 

+ ... 

We define 

Definition 12. yl function K{1) is reproducing for a given lattice A if for any pair of sites i,j we have 

^is:(dist(i,m))i^(dist(TO,j) < Aii'(dist(i, j)), (182) 

m 

for some constant A. 

For a square lattice in D dimensions and a shortest-path metric, a powerlaw K(l) ~ is reproducing for sufhciently 
large a. An exponential decay is not reproducing. However an exponential multiplying a sufficiently fast decaying 



power is. Using this definition and lemma (16 1 and Eq. (181), we arrive at the bound for a reproducing K{1) that 



G{t,j) < 2X(dist(*,j))(2|5| + ^^A+^^A2 + ... (183) 



< 2X(dist(i,j)) 



2! 3! 

exp(2A|.s|) - 1 



A 

Now, we consider a specific case. In order to maintain generality in our choice of F{t), we want one more definition. 
Definition 13. A function f{x) is 0/ decay class g{x) if 

fix) < cig{c2x) (184) 

for all a: > 0, for some constants Ci, C2. 

In some cases we will specify the constants. 

Consider a function g{t) defined to have \g{t)\ < 0{exp{—\t\e{\t\))), as constructed in (TH]. Suppose e{t) is any one 
of the following functions: 

(185) 



\0g{tr ' log(i)log(l0g(t))2 ' I0g(t)l0g(l0g(t))l0g(l0g(l0g(t)))2'- 
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From any such function, in the first se ction of the appendix we succeed in constructing a function F{t) which is of 
decay class g. One may note that Eq. (161) has the extra factor log(2 + \t\)'^ multiplying the decay of F{t), but this 



is still asymptotically bounded by C'(exp(— C2|i|/ log(2 + |i|)^)), for any C2 < 1, so indeed F{t) is of decay class g{t). 

We now restrict to finite dimensional lattices, so that the number of sites within distance I of any site grows at 
most polynomially with I. Then, we find that the function E is also of decay class g. However, here it is important 
to specify the constants: the constant ci is proportional to I/7, while the constant C2 is proportional to j/vlr- The 
function K is of decay class E with constants ci , C2 that depend on the dimension: the number of sites m in the sum 



in Eq. (16) will be dimension-dependent. We are not interested in these dimension dependent constants, but only in 
the constants that depend on 7. Note that K is of decay class g with ci proportional to J'/7 a-^d C2 proportional to 
i/vlr- 

The functions K which arise from these functions g decay faster than 1/a;" for any power a. So, for any a > 0, for 
any C2 < 1, K{x) can be upper bounded by 

K{x) < K'{x)/x", (186) 

where K'(x) = ciK{c2x) for some ci. Note that for sufficiently large a such functions K' are reproducing, and indeed 
one can pick A < C{J' /j), for a dimension-dependent constant C of order unity. Also, K' is of decay class K. Define 

G{l)^ma4f'-''^-'G{i,j), (187) 
so that G{1) is an upper bound on G{i,j) if dist(i, j) = I. Then, 



Lemma 17. For any of the functions g{t) arising from e in Eq. [185), the corresponding G{1) is bounded by 
exp(C|s| J'/7) — 1 times Cig{c2l^ /vlh), for some constant ci,C2 which are dimension dependent, but do not depend 
on 7 or vlr, arid 

Similarly, in any finite dimensional lattice, after we sum over sites i,j we find that 



Lemma 18. For any of the functions g{t) arising from e in Eq. (185), the commutator \\\U sO aUI , O b\\\ is bounded 
by \A\ times by exp(C|s| J'/7) — 1 times Cig{c2l"f /v^ji), for some constants c'i,c'2 which are dimension dependent, but 
do not depend on ^ or vlh. 



Optimization 



In the next section, we will apply this lemma ( 18 ) to specific problems, where the error estimates depend upon 
the error in lemma ( 18 ) for / of order the system size L. Hence, we will find errors that decay with system size L 
similarly to the decay of g{t) with t. In this section, we discuss optimizing these error estimates in a little more detail; 
the reader may prefer to skip this until later. The reader will perhaps realize that instead of choosing a function, 
such a exp(— log^(t)), we could easily have chosen the function exp(— 2|t|/ log^(t)) or indeed exp(— C|f|/ log^(i)) 
for any constant C. One may then start to get optimistic: if, for example, we prove that a certain error decays as 
exp(— lOOOi/ log'^(L)), surely this is much better than proving that it decays only as exp(— L/ log^(L)). Of course, 
the trouble is the constant in front that we have not written: all the of the bounds on the decay of g{t) are up to 
a constant in front, and by choosing a g{t) which decays asymptotically as 0(exp(— 1000|t|/ log^(i))), we will find a 
much worse constant than if we had chosen g{t) to decay only as as C'(exp(— log^(t))). In fact, one will find that 
one gains only for very large t, and hence, for very large L. 

The next possibility one may consider is: suppose, for definiteness, we choose g{t) to decay as exp(— C|t|/ log^(i)) 
for a constant G. How can we optimize the constant C to obtain the best decay for a given system size L7 After 
all, perhaps if L is large one should choose a large constant G in defining g{t), while for smaller L perhaps a smaller 
constant is better. However, as soon as one starts trying to optimize over constants, one should also try to optimize 
over functions. In order to obtain the best possible bounds, one would like to be able to consider all possible g{t) for 
a given L. Thus, the problem which one would really like to solve is: for a given distance L, and a given dimension D, 



what function g{t) will minimize the quantity exp(C|s| J'/7) — 1 times c[g{c2Lj/vLR) in lemma (18) for |s| of order 
unity? In particular, how does the minimum over g{t) of this quantity behave for large L? It should be clear that the 
term exp(C|s| J'/7) — 1 is important. While |s| is of order unity, if G becomes large (for example, L-dependent) then 
this term will become large. 
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Applications 

These operators immediately simplify and tighten several results. For example, |21l I22j rely on quasi-adiabatic 
continuation operators; using the operator here, the superpolynomial decay of the ground state splitting in those 
papers can be tightened to subexponential. 

Further, these operators can be directly inserted into the proof in [S]. This greatly simplifies the construction of 
the quasi-adiabatic evolution operator. Further, one immediately finds that 



Theorem 4. Let g{t) — exp(— te(t)), where e{t) is any function in Eq. (185). For Hamiltonians on a torus, as 
considered in IQj, for fixed R, J/7, the Hall conductance is quantized to an integer up to an error which is bounded by 
a function e{L), where e is of decay class g. 



* Electronic address: mahastin@microsoft . cornj 
[1] P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] E. H. Lieb and D. W. Robinson, Commun. Math. Phys. 28, 251 (1972). 
[3] M. B. Hastings and T. Koma, Commun. Math. Phys. 265, 781 (2006). 
[4] B. Nachtergaele and R. Sims, Commun. Math. Phys. 265, 119 (2006). 

[5] M. B. Hastings, "Locality in Quantum and Markov Dynamics on Lattices and Networks", Phys. Rev. Lett. 93, 140402 
(2004). 

[6] D. M. Basko, I. L. Aleiner, and B. L. Altshuler, "Possible experimental manifestations of the many-body localization", 

Phys. Rev. B 76, 052203 (2007). 
[7] C. K. Burrell and T. J. Osborne, "Bounds on Information Propagation in Disordered Quantum Spin Chains", Phys. Rev. 

Lett. 99, 167201 (2007). 

[8] M. B. Hastings and S. Michalakis, "Quantization of Hall Conductance for Interacting Electrons without Averaging As- 
sumptions", arXiv:0911.4706 
[9] M. B. Hastings, "Lieb-Schultz-Mattis in Higher Dimensions", Phys. Rev. B 69, 104431 (2004). 

[10] M. B. Hastings, "An Area Law for One Dimensional Quantum Systems", JSTAT, P08024 (2007). 

[11] M. Aizenman and G. M. Graf, "Localization Bounds for an Electron Gas", J. Phys. A 31, 6783 (1998). 

[12] M. B. Hastings and Xiao-Gang Wen, "Quasi-adiabatic Continuation of Quantum States: The Stability of Topological 
Ground State Degeneracy and Emergent Gauge Invariance", Phys. Rev. B 72, 045141 (2005). 

[13] M. B. Hastings, "Quasi- Adiabatic Continuation in Gapped Spin and Fermion Systems: Goldstone's Theorem and Flux 
Periodicity", JSTAT, P05010 (2007). 

[14] E. H. Lieb, T. D. Schultz, and D. C. Mattis, Ann. Phys. (N. Y.) 16, 407 (1961). 

[15] M. B. Hastings, "Sufficient Conditions for Topological Order in Insulators", Europhys. Lett. 70, 824 (2005). 

[16] S. Bravyi, M. B. Hastings, and F. Verstraete, "Lieb-Robinson Bounds and the Generation of Correlations and Topological 

Quantum Order", Phys. Rev. Lett. 97, 050401 (2006). 
[17] T. J. Osborne, "Simulating adiabatic evolution of gapped spin systems", Phys. Rev. A, 75, 032321. 
[18] I. Affleck and E. H. Lieb, Lett. Math. Phys 12, 57 (1986). 

[19] A. E. Ingham, "A Note on Fourier Transforms", J. London Math. Soc. 9, 29 (1934). 

[20] J.E. Avron and R. Seiler, Phys. Rev. Lett. 54 259-262 (1985); Q. Niu, D. J. Thouless, and Y.-S. Wu, Phys. Rev. B 31, 
3372 (1985). 

[21] S. Bravyi, M. B. Hastings, and S. Michalakis, "Topological quantum order: stability under local perturbations". 
arXiv:1001.0344, 

[22] S. Bravyi and M. B. Hastings, "A short proof of the stability of topological quantum order under local perturbations". 

[23] J. Dziubanski and E. Hernandez, "Band-limited wavelets with subexponential decay", 41, 398 (1998). 

[24] While this work was in progress, a similar definition was mentioned to me by A. Kitaev, under the name "locally definable" . 



